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obtained for the 2 dimensional 3-state Potts model with a boundary. Exact results are 
obtained on phases with intermediate mobilities. These correspond to non-trivial bound- 
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of Dirichlet and Neumann ones and unlike those trivial boundary conditions, are not 
invariant under a Heisenberg algebra. 
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I. INTRODUCTION 

Conformal field theory (CFT) with boundaries finds applications both to open string theory and to 
various quantum impurity problems in condensed matter physics. These generally describe gapless bulk 
excitations interacting with some localized degrees of freedom. In these problems the gapless bulk exci- 
tations can be represented by a conformal field theory in (1+1) space-time dimensions, often simply free 
bosons. It is generally found that the boundary dynamics renormalize, at low energies, to a conformally 
invariant boundary condition (b.c). 

Quantum Brownian motion (QBM) provides an intriguing example of such a problem. Here one 
considers a heavy particle moving in a d-dimensional periodic potential and interacting with a heat bath. 
A simplified model for the heat bathBij is an infinite set of harmonic oscillators coupkd linearly to the 
particle co-ordinate. (This is the obvious generalization of the Caldeira-Leggett modelD for a particle in 
a double well potential.) 

When the oscillator spectral weight vanishes linearly at low frequencies, corresponding to ohniic dissi- 
pation, the set of oscillators may be represented by a (1+1) dimensional quantum conformal field theory 
of free massless bosons living on a fictitious half-line with the heavy particle at the origin. The number of 
bosonic fields required is given by the dimensionality of the space in which the heavy particle is diffusing 
(normally three). The boson fields at the origin, 0(0), correspond to the momentum of the particle and 

the dual fields, 0(0) to the particle's position. The boson CFT may be regarded as compactified on the 
lattice on which the heavy particle is diffusing. The dimensionless compactification radius (scaled by 
the coupling to the bath) is a crucial parameter for QBM. The fictitious extra dimension is analogous 
to the coordinate along the string, cr, in open string theory, while the field plays the role of the string 
coordinates in D dimensional space time, X^. Conformal invariant boundary conditions obtained in the 
dissipative quantum mechanics framework have therefore immediate applications to open string theory, 
and potential interpretations in terms of D-branes.H 

Until recently it was generally believed that, depending on the strength of the dissipation relative 
to the period of the potential, the particle could only be in either localized or freely diffusing phases 

corresponding to Dirichlet (D):0(O) — constant, or Neumann (N): 0(0) = constant, boundary conditions 



in the CFT, respectively. However, it was recently showro that, for certain lattices, other phases are 
possible with intermediate mobility, in which the particle is neither perfectly localized nor freely diffusing. 
These correspond to non-trivial conformally invariant boundary conditions in the free boson CFT which 
are neither D nor N. In fact, the existence of such phases was shassai earlier in another quantum impurity 
problem: tunneling though a single impurity in a quantum wircHa 

While these interesting phases seem to cry out for a general boundary CFT solution this has, so far, 
eluded us. The problem of classifying all (or all physically relevant) conformally invariant boundary 
conditions in c = 2 conformal field theory (2 free bosons) remains very much open. In the much simpler 
case of 1 free boson, c = 1, it is widely believed that only D and N phases generally occur. On the other 
hand, for conformal field theories with a finite number of conformal towers, an elegant way of generating 
non-trivial conformal boundary conditions is provided by fusion. This method can be used to generate an 
apparsatly complete set of conformally invariant b.c.'s in, for example, the Ising modellJ or 3-state Potts 
modelJlil In the case of Wess-Zumino-Witten (WZW) models, one can obtain by fusion the-boundary 
conditions that describe the low temperature phases of various generalized Kondo problems. l£I However, 
the infinite number of conformal towers present in c = 2 CFT makes the generalization of the fusion 
technique rather subtle. Indeed, it is not even clear in general how to understand the simple D and N 
b.c.'s this way. Whether or not some generalized notion of fusion will ultimately provide a complete 
solution to this problem for c = 2 is unclear at present. Needless to say, the generalization to larger 
values of c, in particular £=3 corresponding to 3-dimensional QBM, is also an open problem. 

The original argumentoO for the existence of non-trivial phases came from a study of the phase dia- 
gram in the quantum wire problem. This was studiedtl by a type of "e-expansion" . By fine-tuning the 
compactification radius the non-trivial fixed points can be moved arbitrarily close to either N or D fixed 
points, allowing a perturbative expansion. The other solved case for c=2 was obtained in the context 
of QBM on a triangular lattice. For a special value of dissipation strength (compactification radius) a 
beautiful exact mapping of the problem onto the critical point of th|e-|3-channel Kondo problemEj was 
obtained by Yi and Kane,Ll allowing use of results from that problenJla based on boundary CFT in the 
SU{2)3 WZW model. Although it can be seen that non-trivial phases exist for a range of dissipation 
strength on the triangular lattice, no general solution has been found. _ . .. . 

Several years after the Caldeira-Leggett type model was proposed for QBMa it was arguedtJO that 
this model does not provide a valid description of a heavy particle interacting with phonons or electron- 
hole paics. In particular, the localized phase does not occur in physical models of quantum Brownian 
motion.Ej Thus physical applications of the model may only be to the quantum wire and related prob- 
lems. Nevertheless, the Caldeira-Leggett type QBM model is perfectly consistent as an ultraviolet (UV) 
regularization of a boundary conformal field theory. 

In this paper we re-examine the soluble triangular lattice QBM Yi-Kane model. Our purpose is both to 
understand this fascinating system better and also to shed light on the general boundary CFT problem. 
Rather than relating this problem to the SU{2)3 WZW model we instead relate it to the 3-state Potts 
model. This is done via a conformal embedding whereby the bulk degrees of freedom are represented by 
a direct sum of the 3-state Potts model, an Ising model and a tri-critical Ising model, satisfying: 

C = Cpotts + Clsing -I- Ctricritical = 4/5 + 1/2 + 7/10 = 2. (1.1) 

The Z3 symmetry of the Potts model corresponds to the point group symmetry of the appropriate 
triangular lattice model. This seems like a more natural formulation of the QBM problem since it does 
not contain an SU(2) symmetry. The 4 conformally invariant boundary conditions in the Potts model: 
free, fixed, mixed and "new" are shown to correspond to the 4 phases that can occur in the QBM problem. 
In this way it is possible to obtain the various fixed points by fusion in the Potts sector of the theory. 
These four fixed points correspond to localized and freely diffusing fixed points, the non-trivial fixed point 
of Yi and Kane and one new fixed point not previously known. The mobilities of these phases and also the 
groundstate entropies can be calculated from the fusion approach. We also show that it is possible to use 



an alternative conformal embedding to Eq. (1.1) where the Ising and tricritical Ising sectors are replaced 
by a Zg conformal field theory, which is a sort of multicritical Potts model.tZl We demonstrate that 
some of the crossovers between these fixed points induced by relevant boundary operators can be studied 
exactly using integrability techniques. This would allow, for example, exact calculations of universal 
quantities at all temperatures for a system which is crossing over between freely diffusing behavior at 



higher temperatures and locaHzed or "non-trivial" behavior at T = 0. The integrabihty also helps to 
understand the phase diagram of the model by determining all RG flows from a given fixed point. 

In the next section we review the fusion approach to boundary CFT and the standard D and N b.c.'s 
for free bosons. We also show there that it is possible to obtain the D b.c. from the N b.c. by fusion 
with a twist field for the c=l case of a single periodic boson, for certain rational radii. In Sec. Ill 
we introduce the QBM problem and its connection with boundary CFT, with careful attention to the 
boundary conditions and the groundstate degeneracy. We show that this degeneracy scales as the size of 
the space on which the particle moves. We derive, apparently for the first time, the connection between 
the particle co-ordinate and the value of the dual field at the origin. In Section IV we introduce a model 
of QBM on a triangular lattice which generalizes the model of Yi and Kane, and conjecture its phase 
diagram, suggesting the analogy with the Potts model. In Section V we discuss our conformal embedding 
and show how fusion in the Potts sector gives the various fixed points conjectured in Sec. IV. We also 
calculate the groundstate entropies and mobilities of these fixed points and discuss the corresponding 
boundary states. We also discuss some additional fixed points, that occur for more general Hamiltonians, 

... (5) 

and are obtained by fusion in the Ising, tricritical Ising or Z\ sectors. In Section VI we discuss the 
integrable flows between the various fixed points. In Section VII we discuss generalizations of this model 
to (n — 1) bosons, corresponding to QBM in {n — 1) dimensions and also discuss the relationship with 
the n-channel Kondo problem. 



II. BOUNDARY CONFORMAL FIELD THEORY 

As will be shown in the next section, QBM can be formulated in terms of (1+1) dimensional massless 
bosons on the half-line, a: > 0, with interactions only at the boundary, a; = 0. As such it is related to a 
number of other problems in particle physics and condensed matter physics that have been successfully 
studied using boundary conformal field theory (BCFT). The basic assumption in this approach is that the 
boundary dynamics renormalize, at low energies, to an effective conformally invariant b.c. A constructive 
approach to classifying possible fixed points is then to attempt to enumerate all possible conformally 
invariant b.c.'s corresponding to a given critical bulk theory. This has met with great success in situations 
where the bulk theory can be formulated in terms of a finite number of conformal towers, with the set 
of conformally invariant b.c.'s being generally in one to one correspondence with the bulk conformal 
towers. In this case, the "fusion" technique is very useful in constructing the conformally invariant b.c.'s. 
Unfortunately, many interesting quantum impurity problems apparently cannot be formulated in terms 
of a finite number of bulk conformal towers, generally corresponding to integer values of c. Consequently 
the problem of studying the fixed points in these cases remains very much open. As we will show in Sec. 
V, for a special choice of parameters, QBM on a triangular lattice can be reduced, using a conformal 
embedding, to a finite number of conformal towers. 

In this section we review two quite different approaches to boundary conformal field theory. The fusion 
approach of Cardy which has been very successful for CFT's with a finite number of conformal towers 
and the straightforward D and N b.c.'s for free bosons. It is not known, in general, how to apply the 
fusion approach for free bosons but we demonstrate that this is possible for a single free boson at a special 
(discrete, infinite) set of compactification radii. 

A. Fusion Approach to Boundary CFT 

We review briefiy some important results due to Cardy on the case of a finite number of conformal 
towers. 

Precisely what is meant by a conformally invariant boundary condition? Without boundaries, confor- 
mal transformations are analytic mappings of the complex plane: 

z = T + ix, (2.1) 

into itself: 



z^w{z). (2.2) 

We may Taylor expand an arbitrary conformal transformation around the origin: 

u;(z)=Va„z", (2.3) 



where the a„'s are arbitrary complex coefhcients. 

They label the various generators of the conformal group. It is the fact that there is an infinite number 
of generators (i.e. coefficients) which makes conformal invariance so powerful in (1+1) dimensions. Now 
suppose that we have a boundary at a; = 0, the real axis. At best, we might hope to have invariance 
under all transformations which leave the boundary fixed. This implies the condition: 

w{t)* = w{t) (2.4) 

where t denotes imaginary time. 

We see that there is still an infinite number of generators, corresponding to the a„'s of Eq. (E|3) except 
that now we must impose the conditions: 

a*^^an- (2.5) 

We have reduced the (still oo) number of generators by a factor of 1/2. The fact that there is still an oo 
number of generators, even in the presence of a boundary, means that this boundary conformal symmetry 
remains extremely powerful. 

Boundary conformal invariance implies that the momentum density operator, T — T vanishes at the 
boundary. This amounts to a type of unitarity condition. Since T(t,x) = T{t + x) and T{t,x) = T{t — x), 
it follows that 

f{t,x)^T{t,-x). (2.6) 

i.e. we may regard T as the analytic continuation of T to the negative axis. Thus instead of working 
with left and right movers on the half-line we may work with left-movers only on the entire line. 

When the bulk theory contains a conserved current, J, it is possible to impose additional symmetry 
requirements on the boundary conditions of the form: 

J(0) = ± J(0). (2.7) 

In the case of free bosons the generic symmetry consists of a produce of U(l) current algebras (one for 
each boson) known as a Heisenberg algebra. It is then possible, more generally, to impose an invariance 
of the form: 

J(0) = 7^J(0), (2.8) 

where the vector J contains the U(l) currents and 7?. is a rotation matrix. Since the chiral energy density 
operators, T and T can be written in Sugawara form: 

T=^P, T^l-f, (2.9) 

it follows that Eq. (|2.6| ) is still obeyed. One of the lessons of this paper is that, while it is fairly 
straightforward to construct boundary conditions in free boson theories that are invariant with respect to 
the full Heisenberg algebra, it is much more difficult, apparently, to find the other boundary conditions 



which obey only the Virasoro condition of Eq. (2.6) and not the additional constraints of Eq. (2.S 



To exploit this symmetry, following Cardy, it is very convenient to consider a conformally invariant 
system defined on a cylinder of circumference P in the r-direction (imaginary time) and length I in the 
X direction, with conformally invariant boundary conditions A and B at the two ends. [See Figure (p.] 
From the quantum mechanical point of view, this corresponds to a finite temperature, T ~ 1/(3. The 
partition function for this system is: 



Z^B=tre-'^^^«, (2.10) 

where we are careful to label the Hamiltonian by the boundary conditions as well as the length of the 
spatial interval, both of which help to determine the spectrum. (We sometimes refer to this as the open 
string channel). Alternatively, we may make a modular transformation, t ^^ x. Now the spatial interval, 
of length, /3, is periodic (the closed string channel). We write the corresponding Hamiltonian as Hp. The 
system propagates for a time interval / between initial and final states A and B. Thus we may equally 
well write: 

ZAB = {A\e-'"'^\B). (2.11) 



Equating these two expressions, Eq. ( 2.10| ) and ( |2.11 ) gives powerful constraints which allow us to 



determine the conformally invariant boundary conditions. 
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FIG. 1. Cylinder of length I, circumference /3 with boundary conditions A and B at the two ends. 



Due to the condition of Eq. (2.6), in the purely left-moving formulation, the energy momentum density, 
T is basically unaware of the boundary condition. Hence, in calculating the spectrum of the system with 
boundary conditions A and B introduced above, we may regard the system as being defined periodically 
on a torus of length 21 with left- movers only. The conformal towers of T are unaffected by the boundary 
conditions. A, B. However, which conformal towers occur does depend on these boundary conditions. 
We introduce the characters of the Virasoro algebra, for the various conformal towers: 

X„(e-'^^/')^^e-^^?(2'), (2.12) 

i 

where Ef{2l) are the energies in the a conformal tower for length 21. i.e.: 

£;r(20 = y<-^, (2.13) 

where the x^'s correspond to the (left) scaling dimensions of the operators in the theory and c is the 
conformal anomaly. 

The spectrum of H'^g can only consist of some combination of these conformal towers, i.e.: 



ZAB=yn'XBXa{e-^^/'), (2.14) 



where the n^^ are some non-negative integers giving the multiplicity with which the various conformal 
towers occur. For minimal conformal field theories, a runs over the finite set of irreducible representations 
of the Virasoro algebra. For more complicated theories like the ones of interest in this paper, the 
corresponding set is infinite, and it is not always clear what kind or representations to expect on the 
right hand side of this equation. We stress that we are interested in boundary conditions restricted by 
conformal invariance only, so characters of generalized symmetry algebras are not relevant here. 

Importantly, only these multiplicities depend on the boundary conditions, not the characters, which are 
a property of the bulk left-moving system. Thus, a specification of all possible multiplicities, n^B amounts 
to a specification of all possible boundary conditions A. The problem of specifying conformally invariant 



boundary conditions has been reduced to determining sets of integers, n'\g. For minimal conformal field 
theories, where the number of conformal towers is finite, only a finite number of integers needs to be 
specified. 

Now let us focus on the boundary states, \A). These must obey the operator condition: 



[T{x)-T{x)]\A) =0 (Vx). 



(2.15) 



(Note that, after the modular transformation, x denotes the periodic co-ordinate.) Imposing periodic 
boundary conditions along the boundary and Fourier transforming with respect to x, this becomes: 

[L„-L_„]|A)=0. (2.16) 

This implies that all boundary states, \A) must be linear combinations of the "Ishibashi states" :t3 

\a) =^\a;m) (g) \a;-m). (2.17) 



Here m labels all states in the a*'' conformal tower. The first and second factors in Eq. ( 2.17 ) refer to 
the left and right-moving sectors of the Hilbert Space. Here, a belongs to the set of representations of 
the Virasoro algebra which appear simultaneously in the right and left sector of the theory with periodic 
boundary conditions. It is in general a smaller set than the one appearing in Eq. (2.11). Thus we may 
write: 



\A)=Y,\a){aO\A). 



(2.18) 



Here, 



|aO) = \a;0)(g)\a;0). 



(2.19) 



(Note that while the states, \a;m) <Si \b;n) form a complete orthonormal set, the Ishibashi states, \a) 
do not have finite norm.) Thus, specification of boundary states is reduced to determining the matrix 
elements, (aO|^). For minimal conformal field theories, there is a finite number of such matrix elements. 
Thus the partition function becomes: 



Z^B=5](A|aO)(aO|B)(a|e-'^^| 



(2.20) 



From the definition of the Ishibashi state, \a) we see that: 



(a\e'^"^a 



E' 



-2i£°J/3) 



(2.21) 



the factor of 2 in the exponent arising from the equal contribution to the energy from T and T. This can 
be written in terms of the characters: 



We are now in a position to equate these two expressions for Zab '■ 

Zab = ^(A|aO)(aO|i3)xa(e-*^'/^) ^Y.^'^BXaie-''^^'). 



(2.22) 



(2.23) 



This equation must.bp true for all values of 1/ p. It is very convenient to use the modular transformation 
of the characters :E3E 



Xa(e 



-ttP/I- 



= Y^Slx,{e^'-''')- 



(2.24) 



Here we refer to 5^ is as the matrix of modular transformations. We thus obtain a set of equations 
relating the multiplicities, n'^g which determine the spectrum for a pair of boundary conditions and the 
matrix elements (aO|^) determining the boundary states: 

J2SSn\s = {A\aO){aO\B). (2.25) 

b 

We refer to these as Cardy's equations; they basically allow a determination of the boundary states and 
spectrum. The problem is to find a set of boundary states A (defined by the coefficients {A\aO)) satisfying 
Cardy's equations, such that, for any pair A, B in this set, the 'n\g are non-negative integer coefficients, 
the identity representation appearing at most once, n\g < 1. An important point in deriving (2.25|) 



is the linear independence of the characters. This property may actually not hold for "non-diagonal" 
theories. In such cases, there will generally be som e othe r quantum number that allows one to distinguish 
the corresponding operators, and still prosie-Eq. (^.25). The study of these equations has given rise to 



considerable theoretical activity recently.EiTES In the case of rational conformal field theories, and with 
the additional key hypothesis of completeness, it has been shown that finding sets of conformally invariant 
boundary conditions is equivalent to finding integer valued representations of the fusion algebra. The 
latter problem is under reasonable control, and as a result, conformal boundary conditions for minimal 
models have for instance been classified in Ref. ( P4 ). 

In the present paper, we are interested in a more difficult problem: in fact, since we look for boundary 
conditions respecting only conformal invariance, the two boson problem has the complexity of an irrational 
theory, even for rational values of the radius. Very little is known about this case, and we rely heavily on 
a not pitirely systematic but highly efficient method callcd-Eusion, which is inspired by Cardy's seminal 
paper ,E£l and was used with success in the Kondo problem.E£l 

Generally, boundary states corresponding to trivial boundary conditions can be found by inspection, 
i.e., given n^^ we can find (a|^). We can then generate new (sometimes non-trivial) boundary states by 
fusion, i.e. given any conformal tower, c, we can obtain a new boundary state \B) and new spectrum n'\g 
from the "fusion rule coefficients", -/V^^. These non- negative integers are defined by the operator product 
expansion (OPE) for (chiral) primary operators, 0^. In general the (OPE) of (j)a with (j)^ contains the 
operator ^c -^ab times. We denote by B the new boundary condition and boundary state, generated by 
fusion from the boundary condition A by fusion with the operator c. The partition function, Zbd, with 
an arbitrary boundary condition D at the other end of the system, is given by the multiplicities: given 
by a fusion with a primary field c has multiplicities given by the fusion rule coefficients A'^^ as 

n'hB = J2^>AD, (2.26) 

b 

hence the name "fusion construction" . Because A^^^ are non-negative integers, the new multiplicities n^j^ 
are also non-negative integers, and thus physical. Alternatively, in the closed string channel, the new 
boundary state |i3 > is defined through 

(aO|i?) = (aO|A)§, (2.27) 

where labels the conformal tower of the ide ntity operator. 

Let us now show that the boundary state (^.27 ) gives the multiplicity ( ^ . 26 ) , i.e. they satisfy Cardy's 



equation( 2.25| ). The right-hand side of Eq. ( 2.25| ) becomes 



{D\aO){aO\B) = {D\aO){aO\A)^. (2.28) 

'-'0 

The left-hand side becomes: 

E ^bnoB - E SSNlnU- (2.29) 

b b,d 



We now use a remarkable identity relatins-the matrix of modular transformations to the fusion rule 
coefficients, known as the Verlinde formulae 



b 



oa Qa 



oa 
"^0 



(2.30) 



This gives: 



T.^>DB = §E'^>^^ = ^{D\aO){aQ\A) = {D\aQ){aO\B), (2.31) 



Qa / V cc UJ\ Qa 

"^0 J "^0 



proving that fusion does indeed give a new solution of Cardy's equations. The multiplicities, n^g are 
given by double fusion: 

n%B^J2^ScNLnU. (2.32) 

b,d 

[Recall that \B) is obtained from \A) by fusion with the primary operator c] It can be checked that the 
Cardy equation with A = _B is then obeyed. 

As mentioned before, the key problem in boundary CFT is the construction of a complete set of bound- 
ary states (andb.c.'s), that is the largest possible set of boundary states, \Ai) such that Za^a- is a physical 
partition function. By physical we mean that the partition function with any two boundary conditions is 
always given by non-negative integer multiplicities n\_j^, . Noting that any linear combination of physical 
boundary states with non-negative integer coefficients, 

J2n,\A,), (2.33) 

i 

also gives a physical partition function in this sense, we see that we must impose an additional condition 
to eliminate such states. The lowest energy state with the same b.c. at both ends of the system is 
independent of the b.c. corresponding to the absolute finite size groundstate, with xt = 0. (This follows 
from making a conformal mapping to the half-plane and using the fact that the identity boundary 
operator exists with any b.c.) We may choose to impose an additional condition that ZA^Ai contains 
the zero energy state exactly once, n^.^. = 1. This eliminates linear combination states. Of course, 
once we have found a complete set, we may always form linear combinations with non-negative integer 
coefficients. Actually, the fusioB-co^truction sometimes gives such liiiear combination states. These can 
also arise from integrable flowgH3~EZl and have physical applicationa^ in some cases. (The existence of 
extra dimension boundary operators is often associated with first order phase transitions.) 

In the case of minimal diagonal theories, copfiorpal boundary conditions are in one to one correspon- 
dence with the set of Virasoro representations Ejcj The corresponding boundary states can be generated 
by fusion starting from the the boundary state |0) defined by its diagonal partition function, n?= = (5°°. 

In more complicated cases, fusion has proven an invaluable tool. One of the key advantages of the 
method is that one can use fusion with fields that are not in the spectrum of the bulk theory. Indeed, as 
already explained above, while the the boundary state (in the closed string channel) must be compatible 
with the bulk spectrum, the conformal towers appearing in the open string channel do not have to be in 
the bulk spectrum. In many cases, it is then possible to use not the fusion algebra of the original bulk 
theory but another one, still compatible with the modular transformation matrix, to generate boundary 
states. For instance, in the case of the three state (^4,1)4) Potts model, the characters appearing in 
the bulk partition function are a subset of all the possible characters at c = |. Using the full matrix of 
modular transformations and fusion algebra for the associated diagonal tetracritical Ising model {A4, A^) 
led to the discovery of the new boundary condition in Ref. ( |ll|) . We will discuss in the next section how 
the fusion method can successfully be applied to the case of a free boson. 

Before leaving this section, we should emphasize that complete sets of boundary states form compatible 
sets of states which may be overlapping. For example, in the Ising model the three states corresponding 
to spin up, spin down and free b.c.'s are a complete mutually compatible set. However, one could also 



consider imposing an up or down b.c. on the dual spins in the Ising model. Such states should also 
be compatible with the free b.c. but not with the fixed b.c.'s on the original (non-dual) spins. This 
non-compatibility is presumably related to the non-local nature of the dual spin variables when expressed 
in terms of the original spins. This implies that one couldn't define a physical partition function with 
the original spins pointing up at one boundary and the dual spins pointing up (or down) at the other. In 
a physical application one typically begins with a particular b.c. and then wants to find all other ones 
which are compatible with it. For instance, we may be interested in adding arbitrary local interactions 
near the boundary and finding all b.c.'s to which the system can renormalize. 

B. Dirichlet Neumann, and Rotated Boundary Conditions For Free Bosons 

We review the simplest conformally invariant boundary conditions, namely Dirichlet and Neumann 
boundary caaditions and some of their generalizations, for multicomponent compactified free bosonic 
field theory.EEl The Lagrangian density of the free boson theory is given by 

C^^idJ)^ (2.34) 

where is a c-component boson field. We identify the boson field as 

$r~.(f+2TTu (2.35) 

where w is any vector in the compactification lattice F, which is a c-dimensional Bravais lattice. 

Let us consider the finite size system of length I, with periodic boundary conditions. The standard 
canonical quantization gives the canonical equal-time commutation relations 

[^^'(t,a;),n'=(t,x')] =i(5J'''^(5(x-x'-mO, (2.36) 

m 

where H — (f>. Together with the equation of motion d^j^d^cj) — 0, we obtain the mode expansion of the 
field $: 

7/ \ r 27rx ^ pt 



n=l 



where u belongs to F and p is the conjugate momentum of (pQ. (h.c. denotes Hermitean conjugate.) 
Because 0o should be identified as (fiQ ^ (j)o + 2ttT, the eigenvalues of p, which we label v, belong to the 
lattice F* which is the dual of F. This is defined by the condition: 

u-veZ, (2.38) 

for all vectors -u g F and v G T* . Each component of the vector a„L is an annihilation operator; the a;^^ 
satisfies the commutation relation [a:^^,a'^„^] = 5^^5nm and all the other commutators involving a^^ 
vanish. A similar definition applies to a„ij, which commutes with a^^^ and oP mL- 

Now we can decompose the field (p into chiral components, namely left-mover (p^ and right-mover 0/j: 




?i.(x+) = I + I + 4(2- +^1-+ + E f^^"^"^"^' + 1--) ' (2-39) 

n— 1 ^^ '^ 

^L ~ 0i? has the mode expansion: 




^>,x)=^ + ^u+^^/+Vf^e-^(-+*)-Hh.c.)-Vf^e^^(--*)+h.c.]. (2.40) 

n— 1 ^ ^ '^ n—1 ^ ^ 



The chiral component of the energy- momentum tensor is given by T{x'*') — (9+0)^. Tlie mode expansion 
of the energy-momentum tensor gives the generators L„i of the Virasoro algebra: 

^(^^) = 7fE^™^""""^'^- (2.41) 

Using the mode expansion (2.39), we obtain 

im = 2 X! ■ "™-',i ■ ^iL :, (2.42) 

where :: denotes the normal ordering and 

{-iy^ttnL {n > 0) 

^{27ru+p} (n = 0) (2.43) 

+tV^aL,^L (n < 0). 
The generators Lm for the other chirality are given by 

-£m = -^ E : a„i-i,B. ■ aiR :, (2.44) 
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where 



{~iy/ndnR (n > 0) 
-^{-27ru+p) (n = 0) (2.45) 

+iV^al„R (n < 0). 

An oscillator vacuum |vac) of this theory, which satisfies a^i//{|vac) = 0, is thus characterized by two 
sets of zero-mode quantum numbers u € T and v €T* , where v is the eigenvalue of the operator p. We 
will denote the oscillator vacuum with these quantum numbers as \{u,v)). 

Now let us consider the possible conformally invariant b.c.'s. As was discussed in Section IIA, the 



conformal invariance of the b.c, Eq. (2.16), implies that these are linear combinations of the Ishibashi 
states, defined in Eq. (2.17). Of course, the primary states must be contained within the Hilbert space of 
the system (with periodic b.c.'s.) A physical boundary state must satisfy Cardy's equation (2.25|), which 



gives a strong constraint on the coefficients of the Ishibashi states. 

The (multicomponent) free boson field theory, which we discuss in the present paper, has an infinite 
number of primary fields. The classification of the boundary states for this case is much more difficult 
than that for CFTs with a finite number of conformal towers discussed in Section IIA. In fact, at present 
we are far from the complete classification. Even for the single-component c — I theory, for which the 
Dirichlet/Neumann b.c. and their generalizations are generally believed to form the complete set, we 
know no complete proof of this. Here we just present some simple boundary states for a multicomponent 
free boson theory. Some examples of more non-trivial boundary states for c = 2 will be constructed using 
a conformal embedding in Section V. 

An oscillator vacuum | {u, v)) is naturally a primary state with respect to the Virasoro algebra. However, 
there are an infinite number of Virasoro primaries which are not oscillator vacua. The most general 
boundary states would include linear combinations of Ishibashi states based on these complicated Virasoro 
primaries. Nevertheless, as we will show in the following, the simplest boundary states can be written in 
terms of oscillator vacua and creation operators of oscillator bosons. 

A sufficient (but not necessary) condition for the free boson theory to satisfy the conformal invari- 



ance (2.16) is given by 

{dnL - d-nR)\B) = 0, (2.46) 

for all integer n. This corresponds to the imp ositi on of invariance under the Heisenberg algebra, as 



discussed in the previous subsection, near Eq. (2.7). For n = 0, it reads u\B) — 0. It means that, to 
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satisfy eq. ( 2.46 ), we ca n use only the states built on the oscillator vacua with u = 0. It can be shown 
that the condition (2.46) is satisfied by the state 



|(0,z;)))=exp(- 



l^'^nL 



^IrMv))- 



(2.47) 



This might be regarded as an Ishibashi state with respect to the U{c) current algebra. It is a linear 
combination of (an i nfini te number of) I shibashi states with respect to the Virasoro algebra. 

While each state ( ^.47 ) satisfies ( ^.46 ) and hence the conforr nal iii variance, it does not satisfy Cardy's 
consistency condition. Let us consider a linear combination of ( p. 47 ) 



[^(•^o)) =5iv^exp(ii/-0o)|(O,tT))), 



(2.48) 



for a constant g^ and constant vector (pQ. The diagonal partition function on the strip is given by 

ZNN{q) = 9n' ( 4-?\ yj^' (2-49) 






where the summation is over the the entire dual lattice T*. Modular transforming to the open string 
channel, it reads 



ZNNiq) = (47r)^/2T/o(r).g^ 



\v{q)) 



^^2.„- 



(2.50) 



where Vo{T) is the volume of the unit cell of the compactification lattice F. We note that the volume of 
the unit cell of the dual lattice is given by: 



VoiT*)^l/Vo{T). 



(2.51) 



The factor g" /ivil))'^ is the character of the U{c) current algebra, and is a superposition of the Virasoro 
characters with non-negative integer coefficients. Choosing 



1 



ON 



{4ny/*^n{r) 



(2.52) 



the state (2. 



satisfies Cardy's condition and thus is a physical boundary state. The coiistant go 

actually represents the (general ly fra ctional) "ground-state degeneracy" due to the boundarycj. From 
the diagonal partition function (2.50) in the open string channel, we can read off the scaling dimensions 
of boundary operators occurring with the boundary condition. Namely, the scaling dimensions of the 
boundary operators are given by 



^N 



2'Ku'' + (non-negative integer) 



(2.53) 



where u is an element of the lattice F. 

The physical meaning of the boundary state \N{(J)q)) turns out to be the Dirichlet boundary condition 
on (j) namely (/) = (/>o at the boundary. This can be checked by the calculation of the partition function 
in the open string channel imposing this boundary condition. This b.c. is equivalent to the Neumann 

boundary condition on the dual field, (i0/dx|o = 0. We label the state with respect to the dual field 
for convenience in the following section. We note that it is not clear whether the N boundary state is 
the only solution of the Cardy's condition even within the linear combinations of the bosonic Ishibashi 
states ( |2.47] ). 

A similar construction of the boundary state is possible, starting from the condition 



{dnL +d^nR)\B) = 0, 



(2.54) 
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instead of ( ^.46| ). Again we are imposing invariance under the Heisenberg algebra, Eq. (2/7). The bosonic 
Ishibashi state is given by 



|(^,0)))=exp(+^4^.4^)|K0)). 

n 

A physical boundary state is given 

\D{$o)) = 9D^eyi-p {iu ■ $f^)\{u,6))), 



(2.55) 



(2.56) 



where the summation is over the entire lattice F and 

The diagonal partition function reads, in the closed string channel. 



^DD 



(q) = 9d' 



1 

viq) 



E' 



and in the open string channel 



ZDoiq) 



1 



E^ 



v^/{2it) 



The scaling dimensions of the boundary operators are 



27r 



+ (non-negative integer), 



(2.57) 



(2.58) 



(2.59) 



(2.60) 



where v is an element of the dual lattice T* . The physical meaning of this boundary state is the Neumann 

boundary condition for the field 0, or equivalently the Dirichlet boundary condition for the dual field 0, 

i.e. (j) = (JJq at the boundary. 
In the simplest case, c = 1, 

r = {nR\n e Z}, r* = {m/R\m G 2} (2.61) 

where R is the compactification radius. In this case the degeneracies are given by: 

g^ = l/l-K^/^Vm], go ^ tt^^^VR. (2.62) 

For the multicomponent (c > 1) case, there is a simple generalization of the above construction of 
N and D boundary states. Since the generators of the Virasoro algebra ( 2.42 ) ar e bil inears in the 
creation/annihilation operator a„, a sufficient condition for the conformal invariance (2.16D is given by 



{dnL - TZd^nR)\B) == 0, 



(2.63) 



where TZ is an orthogonal ex c matrix. This corresponds to another sort of covariance under the Heisenberg 
algebra, Eq. (2.8). This includes the N case (2.4(:) which corresponds to 7?. = 1 (identity matrix) and the 
D case ( 2.54 ) which amounts to 7?. = —1. Other choices of the orthogonal matrix gives different boundary 
states. For example, when 7?. is a diagonal matrix with diagonal elements +1 or —1, it represents a mixture 
of the D and N b.c.'s. When 7?. is a rotation matrix, it gives a b.c. corresponding to the QBM under an 
external magnetic fieldBj We postpone the discussion of such a problem to a later publication. 
The zero modes of the boson are restricted by the condition (2.63) with n = 0, namely 



{2Tru + v) =TZ{-27ru + v). 



(2.64) 
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Its solution is given as a linear combination of basic solutions 

(-U, w) = ^ A„(u„,w„), (2.65) 

n 

where A„'s are integer coefficients. For such (m, v), 2itu + v form a new lattice 2'kT. Only the oscillator 
vacua with these zero modes can contribute to the boundary state. 
A possible boundary state is given as 

|7^)=57^El("'^)))' (2.66) 



where the sum is taken over (2.65) and 

CXD 

\{u,v)}} = exp(^ -4^7^4^) !(«,«)). (2.67) 

n=l 

The diagonal partition function for this state Zfm turns out to be identical to that for Zj^^, if the 
original lattice F is replaced by F. The scaling dimensions of the boundary operators and g-iactois are 
also easily given by this replacement. 

For the multicomponent free b oson, we have thus constructed an infinite variety of b.c.'s by the gener- 



alization of D/N as in eq. ( Z.63 ), for the infinite possible different choices of TZ. The diagonal partition 
function for these b.c.'s has a structure similar to that for D or N b.c. and can be written by a non- 
negative linear combination of bosonic (Heisenberg) characters q^ / {'r]{q)Y ■ We emphasize that we have 
just given some simple examples of physical boundary states for the free boson theory; we have not 
exhausted all the possible conformally invariant b.c.'s. 

In fact, more non-trivial "non-bosonic" boundary states are possible in the multicomponent free boson, 
at least in some cases. In Section V we will construct such highly non-trivial boundary states and will 
demonstrate that they are indeed not bosonic; they cannot be understood by a generalization of D and 
N as described above. 



C. Dirichlet and Neumann Boundary Conditions From Fusion: c=l Case 

The problem of finding out all possible conformal invariant boundary conditions is still an open one 
for all but minimal theories. Even the case of the free boson c = 1 is not fully understood yet. In this 
section, we discuss this problem further, and demonstrate how the twist field (of dimension 1/16) allows 
one to connect D and N boundary conditions in some simple rational cases. This gives some justification 
for the fusion procedure to be implemented in the following sections in the case c — 2. Some of the 
following results are related to the more abstract approach in Ref. ( pl| ). 

We thus consider a single compact free boson, corresponding to the one dimensional limit of the model 
reviewed in the previous section with the compactification lattice {nR\n S Z}. We restrict to the case 
B? = p'/27r, p' an integer, (p' = 1 is the SU{2) symmetric point and p' = 2 is the square of the Ising 
model.) Tiitv-tprus partition function of the periodic boson can be expressed in terms of the generalized 
charactersE3El 

-J oo 

K[^^'^=K[_f=K'^if\ = - J2 q("^'^+^)'/2^^ (2.68) 



as 

M-l 



A=0 



Zp,r^Y.\^i ' (2-69) 
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with M = 2p' . This decpmposition reflects the symmetry of the periodic boson at this particular radius 
under a chiral algcbraE3(traditionally denoted by Ap') which is larger than the Virasoro algebra, and is 
generated by, besides the stress energy tensor T, the current j and a pair of vertex operators of integer 
spin p' . The primary fields of the algebra Ap' are vertex operators V\, A = 0, . . . , M — 1, with fusion 
coefficients 

Nxt. = Sx+t.,. mod N (2.70) 

while charge conjugation takes V\ ~f Vm-x- 

For the same periodic boson, we now consider partition functions in the open string channel. Two 
kinds of conformal boundary conditions are known beforehand: the Dirichlet boundary conditions, where 
the value of the field on the boundary is fixed modulo the compactification lattice, (f) = 4>o + 2TmR, and 
the Neumann boundary conditions, where the value of the dual field on the boundary is fixed modulo the 
dual compactification lattice, (fi = (pi^ + m/ R. These boundary conditions preserve the Virasoro symmetry, 
but do not, in general, preserve the chiral symmetry Ap' . In fact, only the set of M Dirichlet boundary 
conditions with (t>o = 2^, A = 0, . . . , i\/ — f do so, with 

Zd{\/2R)d{p/2R) = K\-ti (2-71) 

These Dirichlet boundary conditions are presumably the only Ap' invariant boundary conditions. Intro- 
ducing the Ishibashi states |A)) associated with the module A of this algebra, the corresponding boundary 
states read 

M-i c:^l M-i 2i7v\fi/M 

1^) = E ^1^)) - E ^f^TT^I^)) (2-72) 

where the matrix of modular transformations is given by 

-^A = J^^'"'"''^''- (2-73) 

The general DD partition functions -^_D(A/2fl)D(/^/2fl) can naturally be obtained by fusion from the 
partition function Zo{o)d(o) using the fusion algebra coefficients ( p. 701 ). 

It is possible to write some of the Neumann partition functions in terms of the ^^'a's. For instance 

A/-1 

ZNio)Nio)= E ^i''^ (2.74) 

A=o,A even 

Similarly, 

M-l 

Zn(o)n(-k/2K) == E "^a (2-75) 

A— 1,A odd 

All the Neumann boundary conditions do however break the symmetry Ap' . This is most clearly seen by 
writing the Dirichlct-Neumann partition function, which is independent of the radius, and well known to 
be 

^^^ = f E'?^'"^'^'^''- (2.76) 

' n 

This partition function cannot be expressed in terms of generalized characters. The corresponding closed 
string partition function cannot either, but of course, it is expressible in terms of the bulk c — \ Virasoro 
characters, as it should be since both D and N are possible conformal boundary conditions for the free 
periodic boson: 
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Znd = 



V2viq) 



^E(-i)>~"^ 



(2.77) 



To relate D and N boundary conditions by fusion, it is necessary to relax the constraint on invariance 
under Ap' . Of course, what one should really do, since one is only interested in general in imposing 
conformal invariance on the boundary conditions, not any higher symmetry, is to simply use the fusion 
algebra of the c = 1 theory for arbitrary radius; this however poses technical problems - in particular, 
continuous modular transformations - which are, for the moment, insurmountable. To proceed, the 
empirical strategy we use is to consider "variants" of the theory, with different, finite, extended algebras, 
for which fusion can be implemented and, hopefully, meaningful boundary conditions for our original free 
boson found. The natural candidate that comes to mind is the Z2 orbifold, with partition function (recall 
i?2 = p'/2n; M = 2p') 



Z, 



p'-i 



(V) 

A 



2 p' 



2Ek 



A=l,3 



^i:/»- + (-i: 



"g" 



:^E/' 



-l)"g"' 



(2.78) 



The operator content of this theory is quite different from that of the periodic free boson. The most 
obvious feature is that, due to the -^ —cj) identification, the multiplicity two of vertex operators has 
now disappeared: we get instead operators with multiplicity one, the first p' — 1 terms, which correspond 
to operators which we call $> , of dimension j-r . The next term in the partition function stands for two 

degenerate operators ^1 of dimension p'/4 (the two energy operators in the Ising square case), the next 
for the twice degenerate twist a'^^'^ and excited twist t*^*) operators of dimensions 1/16,9/16 (the spins 
and the tensor product of spin and energy for the Ising square case). Finally the last two terms are 
respectively the characters of the identity and the current (which we denote 9) , which is primary in the 
orbifold algebra. We will denote the corresponding characters (which appear as moduli squares in the 
foregoing formula) Ki and Kg. 

To proceed, observe that, presumably, the boundary states of the periodic boson and the orbifold boson 
are not mutually compatible in general: for instance, a Dirichlet boundary state for the orbifold |Do(</>o)) 
should be a linear combination of Dirichlet boundary states for the non orbifold boson at values ±(/)o, 
and the normalization will make \Do{4>o)) incompatible with \D{(J)q)). There are thus various families 
of boundary states we can envision constructing. To start - and get quickly to the point - we consider 
the particular cases of \D{(f)Q = 0)), |iV(0o — 0)) for the periodic boson. The corresponding partition 
functions, in the open string channel, can be fully written in terms of the generalized characters of the 
orbifold algebra: 



Z 



'D{o)D(a) — 

N{0)N{0) = 

Zdn = 



Kt + Kb 



Ki + Ke 



K. 



(V) 



K 



K 



(2p') 



(8) 



K. 



(8) 



2-^^p' 



(V) 



(2.79) 



{p' even). Formulas ( 2.79 ) suggest that it is now possible to connect D and N by fusion, which we now 
demonstrate. 

To do so, we need the fusion algebra of the orbifold theory: it can be found for instance in Ref.( |33|) 
with the following results {p' even): 



ad) 


X 9 


= 




r(i) 




a(i)x 


ad) 


= 


/+$« 


' 2-^X even 


$A 


fiWx 


a(2) 


— 


E 


A odd *A 




X$A 


even 


^ 




aW 
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a(i) X $A o. 



•dd = 

$« X $(;) ^ 

<i>(V X $(2) ^ 
p' p' 

$A X $^ = 

$A X $A = 

$2p'-A X $A = 



,(2) 
.(1) 



*A-M + *A+M 
1 + 61 + $2A 

$2A 






(2.80) 



T^-^^ we 



where the fusion of ^x operators holds only for fi ^ X, 2p' — A, and A is defined modulo p' . 

We now start with Zdd and fuse with the ct'"'^' operator. From ct^^^ x I = c/^^and a^^^ x 
get the partition function Zjsm, since the two terms in the last equation of ( 2.79 ) are the characters of 
the twist and excited twist operators respectively. 

For p' even, we fuse again with a'^^'> . We need the basic fusion rule a^^^ x a^^^ = -^ + </*!,' + J2\ even ^^^ 
and also the rule deduced from associativity ct^-'^' x r^^^ ~ O + cjyi^,' +X!a even 't'>''^ *^^ fusion gives therefore 
Zj^fN- To summarize: by fusion with cr*^^^ we thus found Zoo ^^ Z^d -^ Zn^. 

Just to illustrate the difference with p' odd, we remark that in the latter case, while the first fusion is 
the same, the next one has to be with the other twist operator, a^'^\ Using the equivalent of (2.80) in 
this case 



X eve 



X even 



(2.81) 



we recover Znn indeed. This time therefore we have Zjjd -^^(i) Zj^u -^^la) -^atat- 

To express the boundary states, it is better now to carry out a more complete analysis within the 
boundary conditions of the orbifold model. We thus consider the proper orbifold Dirichlet boundary 
states, which are generally defined by (0o 7^ 0, ixR) 



l^o(0o)) 



1 

71 



|i5(0o)> + |^(-0o))] 



(2.82) 



The corresponding Dirichlet partition functions are 



ZD„(\/2R)D„(ti/2R) = -ft^A-p + K\+^,, A 7^ /i, A 7^ 2p' - /i 
= Ki + Kg + K2X, A = ^ 



K. 



2A 



2 P' ' 



X = 2p' - fi 



(2.83) 



The only orbifold Neumann boundary states compatible with the orbifold algebra are the "end point" 
ones: (/)o = 0, ^. We will discuss them later. 

We now need the modular transformation matrix, which reads: 



S 



2(-l)^ 



i-iy 



2 
2 

2(-l)^ 
4cos^^ 






2p' 










Vp' 

-i-iy-f 




^^"j'" 






(2.84) 
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where the conventions are that the matrix maps I /, 6', $p, , $a, cr^* ^t^' '] onto the same quantities but 
with i,j exchanged. 



After modular transformations 
1 



^D„(\/2R.}D^{fi/2R) 



/W 



p'-l 



2X1 + 2xe + 4i(-l)^±'^xi?^') + 4 E -^ ^ -^ ^^i"''' 



(2.85) 



Expressions of the various boundary states in terms of Ishibashi states of the orbifold algebra follows: 



\D,{X/2R)) = 



{2p'y/i 



V2\I) + V2\d) + V2{-iy 



' W«) 



l'^' 



(2), 






(2.86) 



Setting A = gives the periodic Dirichlct boundary state, up to a factor of \/2: 

p'-i 



\Dm 



{2p'y 



i^ |/> + \e) + i= {\4^) + 10^^))) + V2Y: I0A> 



(2.87) 



A=l 



Modular transformation gives the other partition functions in the closed string channel 

1 



^D{o)N{a) 



/W 



/p'Ki - ^p'Ke 



^N{0)N{0) 



/v 



Ki + Ke + 2^K^y^ 



From this, we obtain the Neumann boundary state 

1 



|iV(0)) = 



(2p')i/^ 



p'|/)-V^|0) + ^(|0«)-|0?^> 



(2.88) 



(2.89) 



Observe that, because of the particular radius we have chosen, the Ishibashi states which appear in \N) 
all appear in \D ) also: this is a necessary condition for N to be obtained from D by fusion, as follows 
from Eq. ( 2.27D and of course would not hold in general. 

So far, the boundary states we considered are all in the periodic sector - the Hilbert space in the closed 



string channel is the one of the periodic boson. As discussed in Ref. ( 53) [see also Ref. ( 31)], there are 
also eight boundary states which involve, in addition, the twisted sector. These states are 



where 00 is at the fixed points of the orbifold, ^o = 0, ttR, and similar Neumann like states: 

with 00 = 0, t:/2R. 

The partition functions found in Rcf.( ^4|) can be written as 



(2.90) 



(2.91) 



^D„(0o)±,-D„(0o)± 


^Ki 


^-D„(</.o)±,-Do(0o)t 


= Kg 


(0o)e,£'o(7rfl-0o)e' 





(2.92) 



Similarly for the Neumann sector one has 
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^NMo).N.i.n-^„y - K?''' + K^^''' + ■■■ K^.'J. (2.93) 



Finally, combining Neumann and Dirichlet gives 



^D„(0o)£,JVo(</.o)e' 



Kf\ resp. Kf (2.94) 



where the choice depends on the signs e, e' and the values of ^Oi '?^o- 

Corresponding boundary states can of course be written, giving rise to rather bulky expressions we will 
avoid here. Rather, we would like to stress that the orbif old b oundary conditions are all connected by 
fusion indeed. Starting with the first partition function in (E^92), we get to the second one by fusion with 
the field 6*, thanks to 6* x 6* = /, so _Do(0o)T follows from _Do('/>o)± by fusion with 6. Similarly, starting 



/,(0 



again with the first partition function in (2.92) and fusing with (\) J gives the third partition function 



using that d x <\>'^} = (jyV , so \Do{4>o)vi) follows from \Do{'Kr — 0o)^2) by fusion with 4>i/ . 

The Neumann and Dirichlet sectors are still related through cr*-*', or, also r'-*', using 9 x a^^^ ~ r*-*'. 
Finally, fusion with and 0^, gives the remaining partition functions in the Neumann sector. As for the 
continuous orbifold Dirichlet and Neumann boundary conditions, they are obtained by fusion with the 
operators (j)\. 

One of the lessons illustrated here, is that to obtain all possible conformal boundary conditions for 
the free boson theory, one needs to consider its Z2 orbifold as well. It is not clear in general what good 
"variant" of the bulk theory is needed to explore all possible conformal boundary conditions for a given 
model. In the case c = 2 for example, we shall see that fusion with a field of dimension 1/8 allows one 
to go from D to N. This is a natural generalization of the c = 1 case, the 1/8 field being now a "double" 
twist field for the bosons (t>i.2- We will also see however that fusion with a field of dimension 2/5 is also 
necessary to explore the possible boundary states. The origin of this field in the two boson language is 
still mysterious. Is it some other sort of twist field? 

The discovery of non trivial boundary conditions should also lead to new bulk theories. In the c = 1 
case for instance, consideration of the possible D and N boundary conditions should lead one to the 
discovery of the Z2 orbifold, were this theory not already known. This is also related to the existence 
of a c = 1 theory - the square of the Ising model - which is not a periodic boson. In the c = 2 case, 
the "embedding" we will use to describe partition functions with boundaries leads similarly to a bulk 
theory - the product of Ising, tricritical Ising and Potts - which does not seem to be any simple sort of 
two boson orbifold. Understanding its meaning would be crucial to generalize our work to other values 
of the triangular lattice parameter, a. 

III. QUANTUM BROWNIAN MOTION AND CONFORMAL FIELD THEORY 

We first discuss the one-dimensional case. The study of Brownian motion of a classical particle described 
by the simple Langevin equation 



is of course one of the basic topic in non-equilibrium statistical mechanics. In (3.1), 77 is a phenomeno- 
logical friction coefficient, V an external potential, M the particle mass, and ^ is a random force that 
obeys 

(e(t)e(t')> =2ryr<5(t-f') (3.2) 
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More recently, there has been a great deal of interest in trying to describe the quantum behavior of 



systems for which the classical motion would be described by the Langevin equation (3.1). The most 
tractable approach has been to couple the "particle" (which might actually represent something quite 
different, like the phase in a Josephson junction) to a bath (an environment) with an infinite number of 
degrees of freedom, which provides both the friction and the fluctuating force. The simplest example of 
that approach is the Caldeira Leggett modelEl where the coordinate Q is coupled linearly to an infinite 
set of harmonic oscillators, with the Hamiltonian 






QAfc 

9fe 2 



(3.3) 



Here Q and P are the co-ordinate and momentum of the particle; Viff) is an arbitrary potential, with 
V((^ + a) = Viff). The exact distribution of coupling constants Afc, masses tti^ and frequencies Wfe is 
important ]or the 'properties of the particle only through the weighted density of states 

J(w) = JV^^J(^-c^fc), w>0. (3.4) 

k 

In the limit where the number of oscillators is infinite, J{lu) becomes a continuous function. In the 
following we will rest rict to the so called ohmic case Jito) = rjuj; with this choice, it can be shown that 
the Hamiltonian (^^) when treated classically yields, after elimination of the bath degrees of freedom, 
the Langevin equation ( [3.11 ). 

By a canonical transformation on the oscillator bath: 



Afe Qk mk 2 

ruk ujI ' Afc 



Qk^ 2' Pk^^^kPk, (3.5) 



we may rewrite the Hamiltonian in the form 

J2 



H-^ + ViQ) + lj:[§-+mk^l[,k-Qf}. (3.6) 



where: 



rhk = Xl/mkLol 



^k = ^k 
— \2 



Afe = K/'^k^k = mkUJk (3-7) 

are the masses, frequencies and couplings of the trans form ed oscillators. Henceforth, we drop the tilde 



notation on these quantities. The Hamiltonian of Eq. (li) is invariant under the translation 



Q —> Q + a, qk ^ qk + a. (3.8) 

As the particle moves in the periodic potential it drags a cloud of other "fictitious particles" with it which 
don't feel the periodic potential V but are bound to the real particle by a quadratic interaction. Due to 



the translation symmetry of Eq. (3.8), there is a total momentum, Py, which is conserved modulo 2-k j a: 

PT = P + J2Pk- (3.9) 

Explicitly, a basis of eigenstates of H may be chosen such that any eigenfunction gets multiplied by the 
phase 6*^'=° under this translation. Without loss of generality, the crystal momentum, Pc, may be chosen 
to lie in the first Brillouin zone, |Pc| < tt/o. 

After the canonical transformation of Eq. ( p^ ) the weighted density of states becomes: 

Ji'^) = "^ X! ^k(^kS{uj - ujk), uj > 0. (3.10) 



k 
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The properties of the original particle are unaffected by a rescaling of ojk (independently for each value 
of k) provided that A^ is scaled by the inverse factor such that AfeWfc is held fixed. 

We note that an alternative way of writing the Haniiltonian is in terms of the total momentum Pt and 
transformed oscillator co-ordinates, 



Qk = Ik 



(3.11) 
which define a transformed canonical set of co-ordinates. The Hamiltonian in these co-ordinates becomes: 



^^2v_£b)!,.«„i^ 



p1_ 

rrik 



2 /2 
mkUJkQ k 



(3.12) 



A convenient choice for the oscillator bath is the Fourier modes of a one-dimensional free massless 
boson quantum field. In order to begin with a discrete set of oscillators, we may define the field theory 
with vanishing boundary conditions in a box of size 2Z, so that: 



UJk — k, k ~ T:n/2l, n — 1,2, . 



(3.13) 



Here we have set an arbitrary "velocity of light" to 1. Strictly speaking an ultraviolet cut off must be 
applied, fc < A, but the universal critical behavior will be independent of the cut off. 

Note that there is a zero- frequency oscillator in the limit Z ^ oo. The field couples to the particle at 
X = and the boundaries of the box are at x = ±1. The oscillator co-ordinates and momenta can be 
expressed in terms of the boson creation and annihilation operators as: 



Qk 



i^/k/2{ak - al), pk = (at -|- a[.)/V2fc. 



(3.14) 



We may then define the field and conjugate momentum field on a fictitious one-dimensional space, x, 
completely unrelated to the physical space in which the particle moves, by: 



(x)^ 



J -'^PkSmk{x + I), Il{x) = - J j'^Qk sin k{x + I 



(3.15) 



We note that (pix) corresponds to a "non-compact" boson field, so no term linear in x is allowed in its 
mode expansion. 

We see that, with this choice of oscillator bath, in order to obtain ohmic dissipation, we must choose, 
Afc — 2t]/1. From Eq. (3.7) this implies that the oscillator masses are given by: 



rrik = Afc/w^. = 2ri/lk'^. 

The fictitious particles become infinitel y he avy at zero frequency! 
In this case the Hamiltonian of Eq. (|3.6|) can be written: 



(3.16) 



H 



2M 



ViQ) 



dx 



dx 



(Jl+^27jQS{x)y 



(3.17) 



It is convenient to separate 0(x) into its even and odd components since the odd part decouples from the 
particle: 



(l,,,oix) = [(t>{x)±^{-x)]/V2. 



(3.18) 



The even field obeys a Neumann b.c, d(pe/dx = at a; = 0. Thus we may equivalently consider a field 
defined on the positive axis, < x < I, with a Neumann b.c. at a: — and a Dirichlct b.c. at x ~ I 
coupled to the particle. This Hamiltonian is equivalent to Eq. (3.6) with: 



(ifc^fc, \k = 2rill, k = T:{n + 1/2)11, n = 0,1,2,... 



(3.19) 
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We might expect the translational symmetry to lead to the usual band structure in which the discrete 
set of cigenstates of given crystal momentum depends in some smooth way on Pc. However, when there 
is an oscillator, qq, of infinite mass a dramatic simplification occurs. This follows from observing that 
any eigenfunction of crystal momentum Pc, can be written in the form: 

* = e'P^io^l;, (3.20) 



where ip{Q, qu) is invariant under the translation of Eq. (3^). Since po does not appear in H, it follows 
that the periodic function, -0, obeys exactly the same Schroedinger equation as the full wave-function, 
\1/, independent of Pc- This implies that all eigenstates come in degenerate sets, where the degenerate 
members are labelled by Pc- Effectively we can give the infinite mass "particle", q^, all the momentum 
at no cost in energy. This degeneracy is strictly infinite if Q is allowed to range over the entire real 
line. Alternatively, if we restrict Q to an interval of size Na then we expect the eigenvalues to form 
nearly degenerate sets of size approximately A/", the number of allowed momenta in the Brillouin zone. 
In a localized phase we can identify this degeneracy with the M minima of V at which the particle can 
become localized. However we emphasize that this degeneracy is an exact property of H (at infinite A/") 
and does not depend on which phase the system is in. We have dwelled on this point here because we 
will find that, in the limit where the oscillator frequencies form a dense set, the groundstate degeneracy 
can take the more general form gA/" where 5 is a universal (generally non-integer) number of 0(1) which 
does depends on which phase the system is in. This rather unusual behavior of the model is related to 
its essentially unphysical nature as a real model of QBM, mentioned in the Introduction. We elaborate 
on this point further in subsection C. 

In the field theory representation, the conserved total momentum operator and transformed oscillator 
co-ordinates are: 

Pt = P + V^4>e{0), K{x)=^e{x)+2^Q5{x). (3.21) 

The transformed Hamiltonian (keeping only the even part) is: 



^ = — 2M — +yiQ) + -^j/- 



#^V+H- 
ax 



(3.22) 



Much of the literature on QBM proceeds by integrating out the field, (j)e{x) to obtain a non- 
instantaneous effective action for the particle. However, we introduce a different approach here in keeping 
with previous analyses of quantum impurity problems using boundary CFT methods. In our approach we 
rather absorb the particle co-ordinate and momentum into a b.c. on the field, in the effective low energy 
Hamiltonian. This has the advantage that correlation functions of the field are immediately accessible. 
For the QBM problem this allows us to study the back-reaction of the particle on the oscillator bath. We 
note that this dichotomy of approaches also exists in other types of quantum impurity problems. There 
are two limits in which the QBM Hamiltonian reduces to a boundary sine-Gordon model: the strong and 
weak corrugation limits. We consider each in turn. 

A. Strong Corrugation Limit 

In the "strong corrugation" limit the potential is large compared to the dissipation. More precisely, we 
require the energy bands produced by the potential V , ignoring the dissipation, to have large band gaps 
compared to the energy scale of the dissipation, ry/Af. In this limit we expect that band-mixing effects 
due to the dissipation will be small and the wave- functions for the particle-bath system will be linear 
combinations of wave-functions in a particular band, ipniQ, Pc)- Here Pc is the crystal momentum and 
the wave- functions, ipn are periodic: 

MQ + a;Pc)=MQ;Pc)- (3.23) 

The solutions of the non-dissipative Schroedinger equation are e^^'''^ip{Q; Pc) and the periodic functions 
obey the equation: 
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\2M 






+ V{Q)\^{Q-P,)^Ei;{Q-P,) 



Projecting the Hamiltonian of Eq. (3.22) into a fixed energy band reduces it to 



/■' 1 



dx 



n' 



12 



(3.24) 



(3.25) 



Here tn{Pc) is the dispersion relation of the n*^ energy band for the non-dissipative problem. Note that 
(/>e obeys a N b.c. at x = and a D b.c. &i x = 1. The degeneracy discussed above corresponds to shifting 
(j)e{x) by Pc/y^- However, this is inconsistent with the vanishing b.c. at a; = / so the degeneracy is lifted 
by effects of 0(1/1) as expected. Upon expanding e„ in harmonics: 

e„(Pc) = An + Bn cos Pca + ... (3.26) 

and ignoring the (less relevant) higher harmonics, we obtain the boundary sine-Gordon model: 



H = B cos[aPc - ^/r]a(j)e{G)] + / dx- 



dx 



H 



12 



(3.27) 



This result can be confirmed by considering a (single band) tight-binding model.0 Introducing operators 



Cj which annihilate a particle at site j, the position operator is represented as: 



^E: 



J. 






(3.28) 



The transformation to the total momentum operator then corresponds to a unitary transformation by 
the operator: 



U = e 



-iV^Q0(o) 



Explicitly: 



U^[T\{x) + y^Q8{x)]U = n(.T) 



(3.29) 



(3.30) 



The part of the Hamiltonian containing only the particle can be diagonalized in momentum space. Defin- 
ing: 



=^ = ^E 



Co-e 



~W] 



(3.31) 



the particle part of the Hamiltonian is: 



Hq = ^e(p)cj,Cp. 



(3.32) 



We find: 



U CpU — Cp4^^/^0(o), 



(3.33) 



and hence the Hamiltonian is transformed into Eq. (3.25), using \/2<j){0) = 0e(O). In the case of nearest 
neighbor hopping the dispersion relation is: 



e(P) = -2t cos Pa. 



(3.34) 
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The above derivation generalizes this result slightly to the case of several well-separated bands with 
arbitrary dispersion relation. Of course, for sufficiently weak dissipation and sufficiently low temperature, 
we will only be concerned with the lowest band. 
It is convenient to define the parameter: 



R = 



1 



V^a' 



(3.35) 



so that the Hamiltonian is invariant under the translation; 



However, as mentioned above, the boson 
translation. We return to this point below. 

The boundary interaction has scaling dimension 



> (j)e + 27ri?. (3.36) 

is not compact so we cannot identify the field under this 



rja 

"27 



27ri?2- 



(3.37) 



This can be verified, for example, by using the N b.c. to replace 0"(O) by 20^(0) and using the standard 
result that the scaling dimension of cos P(J)l is /3^/87r. The boundary interaction is relevant for rja^ < 27r. 
For rjo^ > 27r, we may think of the hopping term as being irrelevant corresponding to a fiow to the 
localized fixed point. For 770^ < 2tt we expect a flow to the diffusing fixed point. Note that either 
increasing the distance over which the particle must hop or increasing the coupling to the heat bath 
promotes localization. 



B. Weak Corrugation Limit 



The momentum operator Pt and co-ordinate Q, can be absorbed into the field, (j^e, by transforming 
to a dual cpcillator bath, a transformation that was introduced by Fisher and Zwerger from a different 
viewpoints This is done by first diagonalizing the oscillator part of the Hamiltonian, whose even part is: 



Hn = 



2Af 



1 



dx 



dx 



n- 



(3.38) 



This quadratic Hamiltonian can be diagonalized by solving the classical Euler-Lagrange equations: 

d2 



Lj 4>{x) 



dx'^ 



^S""^' 



(3.39) 



This is just the Schroedinger equation with a (5-function potential. The odd wave-functions are unaffected 
by the potential and the even ones can be expressed in terms of a phase shift, 5{k): 



<j)k{x) ex cos[fc|a::| -f 5{k)\, 



(3.40) 



with frequency w(fc) — |fc|. [Here k runs over the wave-vectors of the even modes, given in Eq. (3.19).] 
Substituting into Eq. ( J3.39| ) we obtain the phase shifts: 



coto(fc) = . 

7? 



For small fc, this becomes: 



m 



TT kM 

2+V 



0{k'). 



(3.41) 



(3.42) 
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The allowed wave-vectors are given by: 

kl + 6{k)=TT{n + l/2), n=l,2,3,. 
For large I {I 3> M /rf) and small k {k -^ v/^'I) we obtain: 

k « nn/l, {n — 1,2,3,.. .) 
The mode expansion for 4>e(x) now takes the form: 



<Pe{x) = Jj^PkCOs[k\. 



6{k)], 



(3.43) 
(3.44) 

(3.45) 



and similarly for IIe{x). [The p^'s are the same operators , up to a minus sign, that occur in Eq. ( 3.15| ) 
for wave-vectors corresponding to odd values of n in Eq. ( |3.13 ).] Noting that, for small fc, 



cos(fc|a;| +5) « -sin[fc|(a;| +A//77)], 



(3.46) 



we see that the effect of the (/>e(0)^ term in Hq at low energies, is to impose, approximately a Dirichlet 
b.c, 06 = on the even field. (More accurately, the low energy field vanishes at x = —M/rj but, in the 
low energy effective theory (j){x) does not vary on such short wavelengths so we ignore this shift.) It is 
now convenient to introduce a transformed even field, (j)' {x) obeying the D. b.c. at x = 0: 



4>'{x) 



-jY^^ink 



X Pk- 



(3.47) 



The even part of Hq may be written in terms of 0' as simply: 

We now must consider the term in the Hamiltonian H of Eq. ( |3.22| ) that couples the field to the particle: 

H,nt = -iVv/M)MO)PT- (3.49) 



This may be expressed in terms of the pk operators of Eq. (3.45) using the explicit expression for the 
phase shift in Eq. ( 3.41 ) as: 



^(0) - -^^^vf ^ 



kpk 



-f ^Tf + [kMf 



We see that, at low energies this is given by: 



^'»>»ff> 



(3.50) 



(3.51) 



but we do no t make this approximation here, proceeding to a more exact analysis. The entire Hamiltonian 
of Eq. ( ^.22 ) can now be written: 



-'An 



kpk 



2?7 



l[rj^ + {kMy 



-Pt 



<lk 



p2 

5 + ^™- 



(3.52) 



In order to complete the square we have used a crucial identity: 
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I 



2V 



_. oo 

„2 /_^ 



^^ if + {-KuM/lf 



l:C0th(VM) = ^ + 0(e-2'''A^). 



(3.53) 



and dropped the exponentially small term. [Note that we are being a bit cavalier and using the small 
k result, k w T:n/l for all n. It ca n be shown that the corrections to this approximation only affect the 
exponentially small term in Eq. (3.53).] It is now convenient to Fourier transform to position space, 
using: 



^.-yf £...„.*. 



dx 



(3.54) 



To this end it is convenient to introduce the function: 



n^) 



I 



Or, 2^ 



coskx 



2?? Y Vrf+JkMf 



(3.55) 



Here the sum is over k — im/l with n = 1,2,3,... f{x) vanishes exponentially at large x. Up to 
corrections which are exponentially small in Irj/M, we may write: 



Here Kq is a modified Bessel function. We may thus rewrite H as: 



H =- dx 



f-«^'-^ 



n'- 



ViQ). 



It is now possible to make a further transformation that absorbs Pt into (j)': 

^"{x)^^'{x)+ f dx'f{x')PT. 

J X 

We may now also absorb Q by going over to the dual field. This is defined, as usual by: 

d(j)" 



(3.56) 



(3.57) 



(3.58) 



dx 
d4>" 
dx 



H" 
-H". 



(3.59) 



These equations only determine (p" (x) up to a constant term. This may be fixed by imposing the canonical 
commutation relations and Euler-Lagrange equations on n"(a;) and 4>"{x). The finite momentum part of 
the dual field is: 



^'(-) = -/fE 



2 Y^ (Ik cos kx 

k 



(3.60) 



The zero mode term in (f)"{x), which we write 0q, must be chosen to commute with all finite momentum 
modes of 11" (x) but have a non-zero commutator with the zero mode of 11" (x): 



U'^^-Pt I f{x)^4 = ^p-. 



iV^ 



(3.61) 



Using the commutator: 
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which can be checked from the mode expansion of these fields, wc find 

^0 =-Vv I dyf{y)~4,'{y) - ^Q. 



(3.62) 



(3.63) 



Note that (t>"{x) obeys N b.c.'s at a; = and I. We see that Q has become (part of) the zero-mode of the 
dual field. The Hamiltonian can be written: 



H^'- 



dx 



U" 



df_ 
dx 



V 



- I dxf{x)ct)"{x) 





(3.64) 



Now consider integrating out the high wave- vector components of 0" (x) to obtain a low energy effective 
Hamiltonian. Since f{x) vanishes exponentially on the length scale M/rj^ it follows that once we have 
reduced the wave- vector cut off to a value much les s tha n -q/M we may approximate cf)" {x) by its value 
at X = inside the integral in the last term of Eq. ( 3.64| ), leading to the simplified expression: 



H ^ - 

2 



dx 



n 



112 



dx 



V 






(3.65) 



The degeneracy mentioned above, corresponds to shifting (f)"{x), or equivalently Q, by a constant na 
where a is the lattice spacing. This degeneracy is lifted at finite I as expected. Now let us consider a 
particular choice for V: 



V{Q) = Vbcos(27rg/a). 
Our Hamiltonian then becomes the standard boundary sine-Gordon model: 



H='- 



dx 



H" 



df_ 
dx 



VqCOS 



2ttR4>"{0) 



(3.66) 



(3.67) 



where we have again introduced the parameter R defined in Eq. ( 3.35| ). Recall that we obtained the 
boundary sine-Gordon model only after integrating out high frequency modes, reducing the cut-off to 
r]/M. This process will, in general, introduce renormalization of the interactions. In order to be able to 
ignore this renormalization, the original dimensionless coupling constant must be small. We may est imat e 
how small it needs to be by considering the renormalization of the boundary interaction in Eq. (3.67). 
This boundary interaction has a scaling dimension of 



A = ^ = 2TrR^. 



(3.68) 
Thus, upon reducing the cut-off from its bare value, A to t^/M ^ the effective coupling is renormalized to: 



A 



A \7^/M 



(3.69) 



Requiring the renormalized dimensionless coupling constant to be small after reducing the cut off to rj/M 
gives the condition: 



Vo < (77/MA)1-"/''"'a. 



(3.70) 
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Thus, this boundary sine-Gordon model is only obtained in the weak corrugation limit. The dimensionless 
parameter, ^Jr\a = 1/ R may take any value, however. 

We note that the weak corrugation formulation of the problem is actually equivalent at Z ^ cx), 
under a dual transformation, to a tight binding model interacting with a heat bath with a Lorentzian 
weighted density of states, J(w). This can be seen by starting from Eq. ( 3.52 ) and making the duality 
transformation: 



Pk 



-Qk/k, qk -^ kpk, Pt 



'Q; Q 



The dual Hamiltonian is: 



H 



91 

2r]l 



V{Pt 



Je 



Qk + 



2r; 



llr]"^ + (kM)' 



iQ 



The last term has the same form as our original heat bath term in Eq. 

LUk = fc, nik 
and hence a weighted density of states: 



7 2 2 

k Pk 
3f) with 



(3.71) 



(3.72) 



1 



Afc 



2r/ 



litf + [kMf 



J{u;)=7T/2j2 



2-qk 



l[Tf + {kM)'' 



-5{uj-k). 



Taking the limit I -^ oo this gives: 



J{u) = 



rjuj 



rf + UJ^M^ 



This is ohmic below a cut off scale i]/M: 



J{.)^-. 
V 



(3.73) 



(3.74) 



(3.75) 



(3.76) 



Note that the dissipation strength parameter, rj has been inverted. We may also drop the first term in 
Eq. (B.72) in the limit ^ — > oo so that the particle part of the Hamiltonian is diagonal in momentum 
space. This is then equivalent to a tight binding model with a dispersion relation: 



eiPr) ^ V{Pt). 



(3.77) 



This gives another way of understanding why a boundary sine-Gordon model is obtained in both the 
weak and strong corrugation cases. We note however, that the P^/2r]l term, which gave 11" (x) a zero 
mode, plays quite an important role in the following sub-section. 

C. Groundstate Degeneracy and Boson Compactiflcation 

An interesting quantity in quantum impurity problems is the zero temperature impurity entropy or 
its exponential, the "groundstate degeneracy", g. More generally, we may define an impurity free energy 
by subtracting off the bulk free energy (the term proportional to /), then taking the length of the bulk 
system to infinity. The zero temperature impurity entropy is defined this way, with the T — > limit taken 
after the infinite length limit. Note that if we take T ^ first, before taking Z ^ oo, the degeneracy is 
necessarily integer valued since the spectrum of the finite size Hamiltonian is discrete. However, taking 
the limits in the opposite order, we are effectively dealing with a continuous spectrum and it is entirely 
possible to obtain a non-integer (even non-rational) "groundstate degeneracy" g, independent of system 
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size, g has been argued4p be universal and to always decrease under renormalization group flow between 
boundary fixed points.E2l Thus it plays a role in boundary critical phenomena analogous to that of the 
conformal anomaly parameter, c, in two dimensional bulk critical phenomena. A new feature occurs 
when we consider g in QBM; it turns out to be proportional to the length of the interval on which the 
particle is allowed to move. We emphasize that we are taking the length, I, of the fictitious line interval 
used to define the oscillator bath to oo. Thus the oscillator spectrum becomes continuous. On the other 
hand, we are considering a long but finite physical line interval, of length Ma with J\f ^ 1, on which the 
particle moves. 

To orient ourselves let us first consider the partition function in the case 77 = where the particle is 
decoupled from the heat bath and the periodic potential V is set to zero. We only consider the even 
oscillator modes, with energies 7r(n + l/2)/l, n = 0,1,2, .. .. (Recall that (f)^ obeys a N b.c. at x = and 
a D b.c. at a; = Z so that the allowed momenta are proportional to the half-integers.) Thus the oscillator 
part of the partition function is: 



00 

n[ 

n=0 



1 _ f.-^in+l/2)/lT 



(3.78) 



Taking the limit IT ^ 1 wc obtain: 



g77ZT/6 

Zosc -^ — ^, (3.79) 

corresponding to an oscillator groundstate degeneracy of 

gosc = I/V2. (3.80) 

This must be multiplied by the partition function of the particle. To make this well-defined, we place the 
particle in a box of size J\fa, with periodic b.c.'s. In this case the allowed momenta are P — 2TTm/J\fa 
giving the full partition function: 



^- E 



exp 



2T:m\ 1 



Ma J 2MT 



g7riT/6 



(3.81) 



If we now take T — > 0, holding Af fixed, the particle has a unique groundstate with P = so it makes no 
contribution to g. We note, however, that in the opposite limit, {Afa)'^MT ^ 1, the partition function is 
proportional to Af since the particle may be treated classically: 



^^iT/6 f°^ dP _p2/,,,^ aAf [MT 



aM ^^^PV2MT^!^r^i±^.lT/6^ (3 82) 



V2 J-00 2^ 2 

Remarkably, this factor of Af will persist in the partition function as we take T — > at fixed J\f, once we 
couple the particle to the oscillators. If we now include the periodic potential V, the partition function 
behaves the same way at temperatures small compared to the bandwidth. This follows since we are only 
concerned with low momentum states in the lowest energy band, whose dispersion may be approximated 
by: 

ei(F) « — — , (3.83) 

iv J 2M* ^ ^ 

for some effective mass, M* . 

Let us now couple the particle to the oscillator bath, rj > 0, but, for the moment ignore the potential 
V. Thus we are considering a freely diffusing particle. From the mode expansion of the dual field, (t>"{x), 



which obeys N b.c.'s at x = 0, Z with the zero mode of Eq. (3.61), or more directly from Eq. (3.52), we 
see that the spectrum is given by: 
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' ^>1 — 1 



I 



(3.84) 



Here the n,„'s are the occupation numbers of the finite momentum oscillator modes. In this case the 
oscillator partition function is given by: 



nti 



-nm/lT]-! £ 



ttIT/6 



(3.85) 



r?i— 1 

Thus the entire partition function, in the limit of large I, for the freely diffusing particle is given by: 

2 






/27rTO\ 1 
^ \j7^) 21'qT 



(3.86) 



Now taking the limit / -^ oo first before taking T — > 0, we may replace the sum over m by a classical 
integral over Py: 






dP 



21T J-oo 2^ 



-P^/21-qT _ 



2M 1^ _ 



TiT/6 



2 V TT 



We have thus obtained a groundstate degeneracy: 



9d^ff 



aJ\f pq 



2 V TT 



(3.87) 



(3.88) 



Of course, the precise value of g depends on precisely how we have defined the oscillator bath, in particular, 
the boundary condition at x = Z on the original field (j^i^)- However, ratios of g at different fixed points 
are expected to be independent of these choices, depending on the oscillators only through the parameter 
77. This is related to the fact that the degeneracy may be regarded as a product of factors for each of the 
two boundaries of the system. Comparing this calculation to the case discussed in the previous paragraph 
we see that the degree of freedom associated with Pt has become infinitely massive in the limit / — > 00 so 
that the associated partition function remains in the classical regime all the way down to T = 0, yielding 
the factor of N. 

Now let us c onsid er the effect of the periodic potential, beginning with the weak corrugation Hamil- 
tonian of Eq. (3.67). This has a scaling dimension of 2-11 /■qa?' . If the potential is irrelevant, rjo? < 2tt, 
we expect that g remains unchanged. On the other hand, if it is relevant, ija^ > 27r, we expect a RG 
flow to a different fixed point. The nature of this fixed point may be deduced by assuming that Vg flows 
to 00, pinning (j)" (0) at its minima. On the other hand, (j)" still obeys N b.c.'s at x = /. Thus its mode 
expansion takes the form: 



(3.89) 



(x) = ^am + J j'Y sin[7r(n + l/2)x/l]pk 



The oscillator part of Z is the same as in Eqs. ( 3.78 ), ( 3.79 ) for the case where the bath is decoupled 
from the particle, and the sum over m simply gives a factor of A/", the number of minima of the potential 
where the particle can get localized. Thus we obtain, in the localized phase: 



AA 
''- ^ T2 



(3.90) 



Despite all the transformations and the RG flow that went into this result the answer is intuitively 
obvious. The degeneracy is simply M times that of the decoupled bath, reflecting the M locations where 
the particle can be localized. Note that for the decoupled bath, the original field ^e obeys N b.c. at x = 
and D at X = L On the other hand, at the localized fixed point the transformed dual field, i^" obeys D 
b.c. at X = and N at x = L 
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It is a useful check on our r esults to derive the low energy spectra using the strong corrugation form 
of the Hamiltonian, Eq. ( 3.25| ). In this case the simple limit is the one in which the dispersion relation 
e„ -^ 0. This corresponds to the limit of zero hopping between sites when the particle is localized. We 
expect to renormalize to this fixed point whenever e„ of Eq. ( |3.26 ) is irrelevant. This has dimension 
7)0? jl'K and so is irrelevant for 7)0} > 27r. It is very important that the hopping term is irrelevant in the 
strong corrugation formulation whenever the potential is relevant in the weak corrugation formulation 
and vice versa. This follows since the scaling dimensions are the inverse of each other. Thus the diffusing 
fixed point is stable when the localized one is unstable and vice versa. This is consistent with our 
assumption that the coupling constant flows to cx3 when it is relevant. Ignoring the hopping term, e„, 
the partition function is just that of the decoupled oscillator bath, Eq. (3.78), multiplied by the result 
of integrating over P^. Note that the crystal momentum, P^ is restricted to lie in the first Brillouin 
zone. The number of momenta in the first zone when the particle lives on a line of length J\fa is J\f, so 
we obtain a partition function of J\fZosc, the same partition function (and finite size spectrum) as we 
obtained using the weak corrugation formulation. Now let us consider the case where the hopping term 
is relevant. If en renormalizes to 00, 4>e(fi) gets pinned at one of the minima of £« which we can take to 
lie at (/)e(0) = {2tt /a){Pc + m/ y/rf) for integer m. At a; = Z, (j)e obeys a simple D b.c, (/)e(0 = 0. Thus the 
mode expansion for (f>e{x) contains a winding mode in this case: 



(t>eix) 



(l — x) Pc + 2'iTm/a 



For fixed Pc, the partition function is: 

gTrlT/6 °o 



J j'^sm{TTnx/l)p„ 



E 



m— — 00 



-{Pc+27T7n/ay/2lTr} 



(3.91) 



(3.92) 



We must further sum Pc over the Brillouin zone. The result of the combined sums is a sum over 
Pr — 2TTm/Afa with m summed over all integers. Thus we obtain the same partition function as in Eq. 



(3.87), which resulted from the weak corrugation formulation and hence the same finite size spectrum. 

We also note that the same ratio of g factors could be obtained using a compact boson field with D 
and N b.c.'s. In this case we would identify: 



and 



(j)eix) ^ (peix) +2nR, 



'{x) ^ ^"(x) + —. 



(3.93) 



(3.94) 



the smallest possible compactification radius consistent with the BSG interactions. In fact, although the 
original field (t>{x) is non-compact, as we remarked above, the trans form e d field , (f)" {x) can be thought of 
as being compact. This follows from the mode expansion of Eqs. ( 3.61 ) ( 3.63| ) onc e we impose periodic 
boundary conditions on Q. Since Q is identified with Q + Afa it follows from Eq. (|3.63) that 



'(^)-0"(^)-^- 



(3.95) 



The effective compactification radius for the dual boson is incr eased by a factor of JV. The allowed values 
of the zero mode of the conjugate momentum given by Eq. ( |3.6l| ) become Ilg = —2'kR/N' the correct 
values for the periodic boson radius R/N . The N b.c. that we have been discussing corresponds to an 
N b.c. on a periodic boson. However the D b.c. is not a simple D b.c. for the periodic boson. The 
boundary sine-Gordon interaction of Eq. (3.67) has A/" inequivalent minima. Thus 0"(O) may take any 
of these M inequivalent values corresponding to a "multi-Dirichlet" b.c. 
The corresponding N boundary state is: 



|iV"(0 = O)>=5]^^|(O,AA2^m/i?)» 



(3.96) 
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with 

9n = {^--'''- (3-97) 

The D boundary state for this compactification radius is: 

\D^{4>^~4>^))=glY.e^^'^-'"'^\{nR/N,m (3.98) 

n 

with 

.1/4. /^ 



9l^^"'\jjj. (3.99) 

The multiple-Dirichlet boundary state, \MD), is obtained by summing (f> over the M lattice points where 
the particle can be localized: 

M 
\MD'') = Y^ \D''{(j) = mR)). (3.100) 



Using: 



J\f 
5]e^'™"/-^-AA^„,A^p, (3.101) 



we thus obtain: 

\MD'')^y/N\D{4) = {))), (3.102) 

proportional to the ordinary D boundary state with radius R. Thus: 

gl^jj ^ VJ7gD = tt^/*VJ7r. (3.103) 



Note that both degeneracy factors are simply increased by a factor of \U\f compared to the ordinary 
compact D and N values. We may use these degeneracy factors for each boundary to calculated the total 
degeneracy associated with the partition functions. These give us: 

9d^ff = [g^^r = ^^, (3.104) 

gioc = glgliD^^- (3.105) 



These are exactly the same values obtained from an explicit calculation in Eqs. (3.8S) and (3.9(]|) 



We note that this infinite groundstate degeneracy, in the limit where the particle is diffusing on an 
infinite line, is related to the unphysical nature of the Caldeira-Leggett type model for QBM. The more 
physical applications of the model, to the quantum wire problem, involve compact bosons from the 
beginning, and do not have this strange behavior. 

In general, the change in Ing due to a change in b.c.'s at one end of the system appears to be insensitive 
to the compactness or non-compactness of the boson. This arises from the fact that the compactness only 
appears when one considers properties depending on both boundaries whereas g is a property of each 
boundary separately. This can be checked very explicitly in the boundary sine-Gordon case, where changes 
in g can be coptputed perturbatively to a very large order, and matched against the thermodynamic Bethe 
ansatz result£3. In the perturbation theory, there is no dependence on the compactification at all orders. 

We note that, upon imposing periodicity of the Hamiltonian under (p ^ (j) + 2ttR oi (f) ^ (p + l/R, the 
correct operator content at the D and N fixed points is given by the finite size spectrum of the ordinary 

31 



compact theory (radius R). Thus in order to deduce the operator content (which is consistent with the 
imposed symmetry) from the finite size spectrum, it is convenient to use compact bosons of radius R. The 
main effect of the non-compactness seems to be to increase the groundstate degeneracy in the partition 
function by a factor of Af, a phenomena whose origin can be traced back to the infinite mass oscillator 
discussed near the beginning of this section. In our discussion of the triangular lattice case we shall use 
bosons compactified on this lattice. Most physical quantities, such as the change in In 5, are not affected 
by this compactifi-cation. 

Note that in the case when the N (diffusing) b.c. is stable, R > 1/V27r, g^ < gjj and vice versa. Thus 
the "g-theorem" is obeyed; the RG flow always serves to decrease g. 

IV. QUANTUM BROWNIAN MOTION ON A TRIANGULAR LATTICE 

In this section we first review the phase diagram of the boundary 3-state Potts model. We then 
introduce the triangular lattice QBM problem and conjecture a phase diagram by analogy with that of 
the Potts model. This is substantiated and studied in more detail in the following sections. 

A. Phase Diagram of the Boundary 3-State Potts Model 

The 3-state Potts model is a natural generalization of the Ising model to a discrete spin variable that 
takes on three equivalent values. (We sometimes refer to these as A, B and C) As such it is naturally 
related to QBM on a triangular lattice, as will become clear in the next subsection. Here we review the 
boundary Potts model phase diagram. The classical Potts model Hamiltonian contains nearest neighbor 
interactions of these spins such that the energy is (-1) when two neighboring spins are in the same state 
and zero otherwise. The critical behavior of the two dimensional Potts model is equivalent to that of a one 
dimensional quantum chain. In addition to the classical Potts interaction this also contains a transverse 
field term which permutes the spins on each site between the three states with, equal amplitudes. A 
complete set of conformal invariant b.c.'s for the quantum Potts chain consistglil'E3 of only four b.c.'s. 
One of these is the fixed b.c. in which the quantum spin at the end of the semi-infinite chain is constrained 
to take a fixed value (either A, B or C). A second b.c. is a mixed b.c. where the quantum spin at the 
end of the chain is hopping back and forth between two of the states (eg. A and B). A third b.c. is 
the free b.c. which results from simply terminating the bulk Hamiltonian without otherwise modifying it 
near the boundary. Finally there is a fourth b.c. whose physical interpretation is not obvious, and which 
we referred to as the "new" b.c. 

In Ref. ( |ll|) a phase diagram was conjectured for the boundary Potts model in which both the 
transverse and longitudinal fields acting on the spin at the end of the chain was varied. As stated above, 
when the boundary transverse field is positive (the same sign as in the bulk at the bulk critical point), 
the system is at the free fixed point. If a boundary longitudinal field is now turned on which favors the 
A state, the system renormalizes to the fixed b.c. On the other hand, if this boundary longitudinal field 
has the opposite sign so as to equally favor B and C then the system renormalizes to the mixed b.c. 
It should be emphasized that the mixed b.c. corresponds to the system dynamically jumping back and 
forth between the B and C states. In particular, it is invariant under the Z2 sub-group of Z^ which 
interchanges B and C . On the other hand, the broken symmetry b.c. whose boundary state is sum of B 
and C fixed boundary states is unstable against an RG flow to the mixed b.c. If there is no longitudinal 
field and the boundary transverse field is negative, then the system renormalizes to the "new" b.c. The 
special case where both types of boundary fields vanish is a type of degenerate boundary condition for 
which the corresponding boundary state is a linear combination of the three different fixed boundary 
states. 

It turns out that all of these RG flows have direct analogies in the problem of QBM on a triangular 
lattice, which thus provides a particular realization of the quantum Potts chain. 
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B. Back to Quantum Brownian Motion 



We consider a nat ura l extension of the ID QBM model discussed in Sec. 2 to two dimensions. The 
Hamiltonian of Eq. (3.6) is extended by introducing a separate set of oscillators, q'^ for each component 
of the particle co-ordinate vector, Q", with an identical set of masses and coupling constants: 



H 



2M 



V{Q) + J2 



Pk , ^kOJk 



2mi, 



qk 



Q^k 



(4.1) 



In the ohmic case we may again represent the oscillators by free massless relativistic boson fields, however 
we now get a separate boson, (j>"'{x) for each component Q"". We emphasize that the bosons still live on 
a fictitious one dimensional line. The derivation of the boundary sine-Gordon model in the strong and 
weak corrugation limits carries over directly to the multi-dimensional case. In the strong corrugation 
limit the Hamiltonian is: 



Hs^ 



1 



dx 



n" 



dx 



e[Pc-VvHO)]. 



(4.2) 



We have dropped extraneous ' and e notation. Here e{P) is the dispersion relation (in the lowest band) 
for the dissipationless problem. 0(a;) obeys an N b.c. at x = and a D b.c. at x = I. The c-numbers 
Pc can take any values in the first Brillouin zone. Similarly in the weak corrugation formulation this 
Hamiltonian is: 



^w — 77 



ft^l 



V 



-m 



(4.3) 



(We have again dropped the cumbersome " notation.) 0(a;) obeys N b.c.'s at x = 0, I. 

The observations about boson compactness made in Sec. HI carry over to the two dimensional case. 
We henceforth set: 



V 



1, 



by a rescaling of the lattice spacing. We label the physical lattice on which the particle moves T* 
convenient to impose periodic b.c.'s on the particle co-ordinate, Q, so that we identify: 



(4.4) 
It is 

(4.5) 



where a is any vector in T* and A/" is a large integer. It follows from the derivation in Sec. HI of the 

boundary sine-Gordon model for QBM that (f> may be regarded as compactified on this "coarse lattice" 

with spacing bigger by the factor of Af. Equivalently, we may regard (j) as being compactified on the 
lattice r* with an extra factor oi Af appearing in the degeneracy, g. 

We wish to consider a model of QBM in a periodic potential with hexagonal symmetry. The tight- 
binding model that we consider is a simple generalization of the one introduced in Ref. ( R). We consider 
a triangular lattice generated by the vectors: 



ai = (a, 0), as = (a/2, V3a/2). 



(4.6) 



This model exhibits interesting behavior for any value of the dimensionless parameter a but in this paper 
we focus primarily on the value: 



= 4^/3, 



(4.7) 
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since, only then, can it be mapped into the Potts model. The simplest form of the tight-binding Hamil- 
tonian, before adding dissipation, is: 



H^-tj:^ 



(4.8) 



(^J> 



where the sum is over nearest neighbors. 

In order to understand the connection with the 3-state Potts model it is useful to focus on the symmetry 
transformations that leave fixed the center of one of the triangles on the lattice. These consist of three- fold 
rotations and reflections, referred to as 3m in the international crystal nomenclature. There is actually a 
larger point group symmetry, 6mm, when one considers transformations that hold fixed a lattice point. 
However, this is not relevant for our purposes. It is convenient to decompose the lattice into three 
sublattices. A, B and C such that the nearest neighbours of a point in one sub-lattice are in the other 
two, as drawn in Fig. (0) . Note that a 27r/3 rotation of the lattice maps each point on the A sublattice 
into a point on the B sublattice and similarly B ^ C and C -^ A. Simlarly a reflection about a line 
passing through A points interchanges all B points with C points. Hence we may think of this point 
group symmetry as the permulation group on three objects, S^. This is the symmetry of the 3-state 
Potts model. We will also consider an on-site potential which assigns different energies to the A, B and 
C sub-lattices, va, vb and vc, thus breaking the S^ symmetry. 




FIG. 2. Triangular lattice with A,B and C sublattices marked. Note that the A points also form a triangular 
lattice with spacing \/Sa and orientation rotated by 90*^. The B and C points form a honeycomb lattice. 



To map out the phase diagram let us begin with the strong corrugation formulation in the case t > 
and all Vi = 0. t is relevant for this lattice spacing, (of dimension x — 2/3) corresponding to the boundary 
sine-Gordon Hamiltonian: 



where 



Ho= - dx 



Hint = —t 



COS V 277 



Hs = Ho + Hi, 



\ dx 



COS 



(4.9) 





(4.10) 
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(The fields are all evaluated at a; = but we suppress this argument. We have set Pc to 0.) We expect t 
to induce an RG flow to the perfect mobility phase, corresponding to a D boundary condition, (/)(0) ~ 0. 
This fixed point is stable as we can see by considering the lowest dimension potential term with the full 
symmetry of the lattice, in the weak corrugation formulation, which is: 



Hw — Hf) + Hi, 



(4.11) 



where 



B,.-J^_.l. 



H,nt = -V 



cos 




cos VTI" 



02 - V301 j 



(4.12) 



This has dimension 2 and is irrelevant. Now consider turning on a potential va < 0, which favors the A 
sub-lattice. We do this in the most symmetric possible way, preserving a Z^ symmetry of rotation about 
a lattice point, as well as a mirror symmetry about a lattice link, the 3m subgroup of the original 6mm 
point group. 

Choosing the origin to lie on the A sub-lattice, this potential is: 



H -f H -VA 



cos 



^f^k 



cos 




(4.13) 



This operator has dimension 2/3 at the perfect mobility fixed point, which corresponds to a D b.c. on 

(j) or equivalently, an N b.c. on (p; it is relevant. Thus, for va > 0, we expect a flow to a localized fixed 

point, corresponding to a D b.c. on (f>. The stability of this fixed point can be checked by observing 
that, if VA fiows to infinity, then only hopping between A sub-lattice points is possible. The particle gets 
localized on one of the A sub-lattice sites. This is a stable fixed point since the intra-sublattice is now 
V3a and hence the hopping term gives: 



Hi, 



-t 



cos 



47r02 + cos ^/tT f 02 + \/301 j + cos ^/tT f 02 - V301 j 



(4.14) 



of dimension 2 which is irrelevant. Note that the dimension of the nearest neighbor hopping term is 2/3, 
smaller by a factor of 3 due to the reduced distance of the hop and hence relevant. The effect of the 
symmetry breaking term, va is to stabilize a localized phase on one of the sublatticcs. 

The analogy with the Potts model is quite transparent. The perfect mobility fixed point corresponds 
to the free b.c. in the Potts model. The potential va corresponds to a longitudinal boundary field which 
produces a flow to the fixed (A) b.c. Now consider the case va < 0. If va -^ — oo, the particle is 
localized on one of the B or C sub-lattice sites. Together, these two sub-lattices define a honeycomb 
lattice. However, for finite negative va, this is not a stable fixed point. This can be seen by considering 
hopping on this honeycomb lattice, since the nearest neighbor distance is again a, this hopping term has 
dimension 2/3 and is relevant. Thus, following the logic of Yi and Kane, there must be an intermediate 
mobility fixed point. In the Potts model we add a boundary field which favors either the B or C state. 
This leads to a flow to the mixed b.c. of the boundary Potts model. Thus wc see that the new fixed point 
found by Yi and Kane corresponds physically to the mixed boundary fixed point in the Potts model. In 
the Potts model we think of the boundary spin as fluctuating back and forth between B and C states. In 
the QBM problem we think of the particle as hopping back and forth between B and C sublattices. As 
can be shown explicitly, this state has an intermediate mobility. We refer to this at the "Y" fixed point, 
after Yi-Kanc. 

In fact, we can discover yet another intermediate fixed point in the QBM problem by pursuing this 
analogy further. We now set the symmetry breaking potentials Vi, to but consider a negative hopping 
term, t < 0. This would seem to correspond to a negative transverse boundary field in the Potts model. 
This was argued to lead to a fiow to another fixed point, unimaginatively referred to as the "new" fixed 
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point in Ref. ( ^^)- Thus we expect a new fixed point to occur in the QBM problem with t < 0. We 
refer to this as the "W" fixed point. Such a new fixed point is possible due to the lack of particle- hole 
symmetry of this model. For i > the lowest energy single-particle state has crystal momentum p ^ Q. 
However, for i < the lowest energy single-particle states occur at two inequivalent crystal momenta: 
p = ±(27r/3a, n/y/Sa). These two lowest energy states for the particle are dual to the B and C sub-lattices 
in the mixed phase. In fact, duality in QBM is easily understood, simply corresponding to: 



(4.15) 



Comparing Eqs. ( 4.10 ) and ( 4.13 ) we see that the flow from localized to perfect mobility is dual to the 
flow from perfect mobility to localized on the A sub-lattice, induced by a positive t or va respectively. 
Similarly, the flow from localized to mixed is dual to the flow from perfect mobility to new, induced by 
a negative t or va respectively. 

If both transverse and longitudinal boundary fields are set to in the Potts model we get a phase 
referred to as ^ -I- i? -I- C signifying that the boundary spin remains fixed in any of its 3 possible states. 
The corresponding boundary state is a sum of three boundary states corresponding to the 3 possible 
fixed b.c.'s. Turning on a positive transverse field produces a flow to the free b.c. whereas a negative 
transverse field produces a fiow to the new b.c. The analogue in QBM is to start with t = 0, corresponding 
to the particle being localized at any lattice site {A, B or C). We may, if we wish, include an irrelevant 
intra-sublattice hopping, of dimension 2. Now turning on t induces a flow to either perfect mobility or 
new flxed points. This phase diagram is summarized schematically in Fig. (0). Note that, since only 
localized behavior is possible on the line i = 0, two copies of the Y flxed point must occur, as drawn. 
Also note that, although the flow away from the "localized on A, B or C" flxed point depends on the 
sign of i, that from the "localized on B or C" fixed point does not. This is consistent with the fact that 
the sign of the hopping term cannot be changed by a redefinition of the sign of the electron operators for 
a triangular lattice but can be for a honeycomb lattice. 



localize 
on B or C 





diffusing 
(Neumann) 



localized 
won A, B or C 




localized on A 
(Dirichlet) 



new (W) 



FIG. 3. Phase diagram and RG flows of the triangular lattice QBM model with hopping strength t and potential 

VA- Neumann and Dirichlet b.c.'s are imposed on the dual fields, 4>. As pointed out in Ref. ( [7|) and discussed in 
Sec. VII, the "localized on B or C" phase corresponds to th weak coupling fixed point in the 3-channel Kondo 
problem. 



We a lso co mment briefly on the model with general values of a. The strong corrugation Hamiltonian 
of Eq. (4.10) becomes: 
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Hi, 



cos acpi + cos a 




\/302 



and the symmetry breaking weak corrugation term in Eq. (4.13) becomes: 



H ^ H 



VA 



47r0i 
3a 



2n(i 



3a 



27r(/)2 



cos 



27r0i 
3a 



27r(; 



V3a 



(4.16) 



(4.17) 



The scahng dimension of the hopping term at the (symmetric) locahzed fixed point and of the symmetry 
breaking potential at the N fixed point are: 



An 



a^/2TT 
87r/9a2. 



(4.18) 
(4.19) 



Thus the hopping term is relevant for a^ < 27r and va is relevant for a^ > 8tt/9 indicating that the perfect 
diffusion fixed point is unstable in this region. For va > we again expect a flow to the fixed point where 
the particle gets localized on the A sublattice. The scaling dimension of the intra-sublattice hopping 
term is 3a^/27r so this localized phase is stable in this region and we expect an RG flow from diffusing to 
"localized on A" . On the other hand, ii va < we should consider the stability of the "localized on B or 
C" phase. The scaling dimension of the inter-sublattice hopping term (between the B and C sub-lattices) 
is a^/27r, so neither diffusing nor localized phases are stable over the range: 



87r/9 < a^ < 27r. 



(4.20) 



In this entire range we expect a non-trivial fixed point. This can be studied using an e expansion for 
a^ = 87r/9 -I- e or a^ = 27r — e for < e <C 1. It can be solved exactly at a^ = 47r/3 as we show in the next 
section. In general, however, it remains an unsolved problem. 

It is also instructive to calculate the degeneracy (g) factors at the diffusing {va — 0) and "localized on 
A" phases. We may do this using the results in Sec. IIB for compact bosons, taking into account the 
connection of the QBM problem with compact bosons explained in Sec. Ill and earlier in this section. 
We assume that the space on which the particle moves is periodic with 



Q ^ Q+Afdi, 



(4.21) 



where the primitive vectors a^ are given by Eq. (!.(:). The unit cell area for the fine lattice with spacing 
a is: 



Vo(r*) = \/3aV2. 



(4.22) 



(Note that this is twice the area of a triangle on the lattice; each triangle contains 1/2 point or three 
points, each of which is shared with six triangles.) Thus the degeneracy factor for the D fixed point where 
the particle is localized on any oi A B or C sublattices is: 



9?,D 



7V^/2^ 



(4.23) 



This is the result for the compactified boson with dual lattice F* from Eq. ( 2.57 ) multiplied by the factor 
of A/". The degeneracy for the state where the particle is localized on the A sub-lattice only is: 



go = 5313/3, 

and that for the state where the particle is localized on the B ot C sub-lattices is: 

goD = (2/3)g3£). 



(4.24) 



(4.25) 



These factors of 1/3 and 2/3 can be understood as simply corresponding to the number of sites on which 
the particle can be localized. Alternatively, we may think of the D fixed point as being the same as 3D 
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except that the lattice constant a is increased by a factor of -^/S and the factor M, defined in Eq. ( |4.5| ), 
is decreased by a factor of \fi in order that the area on which the particle moves remains fixed. The 
degeneracy factor for the N fixed point is: 



This is obtained from Eq. ( |2.52| ) for a compact boson, multiplied by the factor of M. Note that the total 
degeneracy for the phase where the particle is localized on ^, B or C sub-lattices is: 

9ioc = g3D9N = ^- (4.27) 

Following the discussion in Sec. Ill for the one-dimensional case, we see that this is simply the number 
of lattice sites times the degeneracy for the decoupled oscillator bath. (The factor of I/a/S that occurs 
in the one-dimensional case gets squared in two dimensions.) We also see that the ratio between fully 
localized and freely diffusing fixed points is: 

^ = 4^. (4.28) 

Thus the flow from fully localized to freely diffusing is consistent with the g-theorem only for: 

On the other hand the ratio of g-factors between the "localized on A" and freely diffusing fixed points is: 

9d 47r 



9N 3%/3a2 ' 
so the flow between freely diffusing and localized on A is consistent with the g-theorem for: 



(4.30) 
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It is instructive to consider a more general model with even less symmetry than that of Eq. (4.12) 



The model of Eq. (4.13) in which we have changed the energy on the A sites, relative to B and C sites 
still has a 3m point group symmetry when we consider transformations holding fixed a lattice point (A, 
B or C). The less symmetric mode l that we now consider has no rotational symmetry whatsoever. This 



model is obtained from Eq. (4.13) (weak corrugation formulation) by letting the coefficients of the 3 
terms be different: 



^V2^W-^i 





Hw = Hq-vi cos V27rW -^i - V2 cos V2tt \ j.i'i + \ 7;<p2 - vs cos V2tt \ j.4'i ~^ \ 7^4'2 ■ (4.32) 



The behavior of this model becomes obvious in the case vi ^ Q , V2 — vj, = Q . In this case the (relevant) 
potential only depends on the co-ordinate Qi. Thus we expect the motion of the particle to be localized 
in the 1-direction but freely diffusing in the 2-direction. Now there is only 1 minimum per unit cell for 
either sign of vi. We may think of the particle as diffusing freely along vertical lines on the lattice. To 
check the stability of this fixed point we should consider the hopping process between vertical lines. Since 
the horizontal distance between points on the A sub-lattice is 1/2 times the lattice spacing (of the A 
sub-lattice) we expec t this hopping term to have dimension 1/2 [smaller by 1/4 than the dimension 2 



hopping in Eq. (4.14).] Thus this fixed point is unstable. 
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V. CONSTRUCTION OF BOUNDARY STATES VIA A CONFORMAL EMBEDDING 

A. Conformal embedding 

Let us consider the QBM on the triangular lattice discussed in Section IV, considering the weak 



corrugation Hamiltonian of Eq. (4.17). Because the symmetry among the sublattices A,B and C is 
broken by the applied potential, the generic translation symmetry of the model is the invariance under 
translations which do not interchange the sublattices. We can impose this symmetry on the system, so 
that we can use the compactified formulation. The compactification lattice F* is thus a triangular lattice 
of nearest neighbor distance v3a- Namely, we introduce the compactification 

$^$+V\ (5.1) 

where T* is a triangular lattice generated by 

(0,\/3a), (3a/2,V3a/2). (5.2) 

Its dual r is generated by (2/3a,0) and (l/3a, l/-\/3a). The allowed vertex operators made from cj) 
are of the form exp(jw- (/>), where v is an element of the lattice T*. As mentioned in Sec. IIIC, the 
operator content upon compactification corresponds to the possible boundary perturbations given the 



imposed symmetry. Note that the Hamiltonian of Eq. (4.17) is the most general one, up to less relevant 
operators, respecting the 3m point group symmetry. 

For va — Q, the corresponding boundary condition is Neumann on 0, while it is Dirichlet on cf) for va — > 
oo. However, there are other nontrivial boundary conditions which are neither Dirichlet or Neumann, as 
discussed in Sec. IV. They correspond to nontrivial phases in QBM. Our problem, then, is to construct 
and to classify the possible boundary states of the c = 2 CFT. Possible boundary conditions may be 
restricted by higher symmetries, such as current conservation at the boundary. However, in the case 
of QBM, apparently the only requirement (at the RG fixed point) is the conformal invariance on the 
boundary and Cardy's consistency conditions. Even if there are extra symmetries in the bulk, they are 
not necessarily respected at the boundary. This allows some extra boundary conditions which would be 
forbidden if the higher symmetry is imposed on the boundary. 

Unfortunately, at present we have no systematic understanding of the boundary states in an irra- 
tional CFT, which has an infinite number of primary fields. While we do not know the solution to this 
fundamental problem, in this paper we analyze some nontrivial boundary conditions for the special com- 
pactification radius a^ = 47r/3. In this case, the irrational c — 2 CFT admits a conformal embedding 
in terms of rational CFTs. For such a rational CFT, we can invoke the fusion construction of boundary 
states, to find non-trivial boundary states for the c = 2 CFTg 

This value of a^ is actually the point where Yi and KaneQ obtained the exact value of the mobility 
at a nontrivial fixed point by a slightly different approach. Our approach leads to a somewhat more 
systematic description of the boundary states, and several additional results including novel boundary 
states. 

The present c = 2 boundary CFT admits a conformal embedding c = 2 = 1/2 + 7/10 + 4/5 = (Ising) + 
(Tricritical Ising) + (Potts). Namely, the partition functions on the strip can be expressed (in the open 
string channel), as 

Zooiq) = [xlxl + x{/2xI/2){Xq + xf) 

+ (Xoxiv5 + x[/2xl/w){X2/5 + Xvs)' (5-3) 

ZNN{q) = iXaxl + x(/2xiv2)(x^ + xf + 2x^/3) 

+ (XoX^/5 + Xv2Xvio)(2xf/i5 + xf/5 + xf/5), (5-4) 

where Xh ' i^ ^^^ Virasoro character of the weight h for Ising, Tricritical Ising and Potts model, 
respectively. (The argument q is omitted.) The two amplitudes have the same factors for Ising and 
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Tricritical Ising sector; the only difference is in the Potts sector. We emphasize that these are the 



partition functions for compact bosons compactificd on the lattice, T* of Eq. (5.2). They do not quite 
correspond to the partition functions occurring in the QBM problem. Nonetheless, as explained in Sec 
IIIC, these are the partition functions which are related to the boundary operator content. We re mark 



that, since the and 3 and also 2/5 and 7/5 partition functions only occur added together in Eqs. (5.3) 
and (|]j), Zdd and Z^n can be expressed in terms of the W-characters of the 3-state Potts model, 
invariant under W-symmetry. 

In fact, we have no proof of these identities. However, we have verified, using MATHEMATICA, that 
the first several (40 to 100) terms in the series expansion agree exactly. In the following, we assume 
this conformal embedding. Because of several additional pieces of evidences, we believe it is indeed 
correct. Namely, as we will show later, we can construct several reasonable boundary states based on the 
embedding. 

We also note that, the torus partition function of the c — 2 theory, even at this particular compact- 
ification, apparently can not be written in terms of Ising, tricritical Ising and Potts models. We are in 
a somewhat strange situation that the boundary CFT admits the conformal embedding while the bulk 
CFT does not. Although this is difficult to understand, it does not pose a real problem in our construc- 
tion. The fusion construction is a m ultipl ication by constant factors of each Ishibashi state appearing 



in the reference boundary state, Eq. ( 2.27 ). These constant factors are given by the matrix elements of 
the matrix of modular transformations. Thus, by construction, the obtained boundary states is a linear 
combination of the Ishibashi states. They are well contained in the Hilbert space of the theory, as long 
as the original reference state is a physical boundary state. 

In fact, a similar phenomenon already appeared for the Dirichlet and Neumann boundary conditions in 
the simpler c = 1 free boson theory, as we have discussed in Section II. There, they are connected by the 
fusion with a twist operator, which is contained in the Z2 orbifold theory. While the free boson theory is 
not equivalent to the Z2 orbifold, the latter appears in the discussion of the boundary conditions in the 
former. 

It is a valuable check to see that the dimensions of boundary operators at D and N fixed points are 



reproduced by the conformal embedding of Eqs. ( |5.3| ) and (5.4). We see that at the D fixed point we have 
only integer dimension operators. (The A = 1 operators presumably correspond to the exactly marginal 

operators, dcj), which are not allowed by symmetry.) On the other hand, at the N fixed points we see that 
there are 6 relevant operators with A = 2/3. These must correspond to the 3 cosine operators in Eq. 



(4.32) together with the corresponding sines. It is interesting to note that 2 of these operators are purely 



in the Potts sector while the other 4 are products of operators in 2 or more sectors. We expect that the 2 



pure operators are the interaction term in the Hamiltonian of Eq. (4.17) and the same operator obtained 
by replacing all cosines by sines, as argued below. 

In fact, there is an alternative conformal embedding, also involving the Potts model, in which the Ising 
and tri-critical Ising models are replaced by the Z^ ' conformal field theoryO (A detailed review of Z^^' 
theory is given in Ref. ^) This may be regarded as a sort of tri-critical 3-state Potts model. Again there 
is a W-algebra and only W characters occ ur i n the conformal embedding. This alternative conformal 



embedding is obtained from Eq. (|5.3D and (5^) by replacing the Ising and tricritical Ising characters by 
the .Z3 characters. Here we use the (conjectured) identity 

XoX^ + x(/2Xr/2 = Xo + 2X2, (5-5) 

XoxJ/s + x(/2Xr/io = 2x^/5 + xl/5- (5-6) 

Here Xa refers to characters in the Z| theory. We note that, again, this identity is only verified up to 
finite order with MATHEMATICA. The equivalence between the Z3 ' model and the product of Ising 
and tri-critical Ising models is only a partial equivalence. Certain sums of conformal characters in the 
two theories are presumably equal but there is not a complete equivalence of all characters. 

Fu rther insight into the conformal embeddings can be obtained by rewriting the interaction term in 



Eq. (4.17) in a more symmetric way, first by reexpressing everything in terms of chiral components, using 
'pL = </*/£ at the fixed point, and by introducing the new chiral fields ^j: 
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with propagators 



$1 = 

$2 = 
$3 = 






H (01L + %/302L 
t (^1^ - V302L 



(5.7) 



The 



($,(z)<i>,H) = -M|ln(z-w) 
($,(z)$,H)^ |z^ln(z_z«), z^j. 

(weak corrugation) interaction term in Eq. ( |4.17| ) now reads, recahing a^ = 47r/3, 

Hint =-Wa(c0S$i+C0S$2+C0S$3) = - (wa/2) (*i + ^J ) 



(5.8) 



(5.9) 



where: 



fc=i 



(5.10) 



^1 corresponds to the fundamental Z3 parafermion field of the Potts model.HEj (This generalizes the 
representation of a Majorana fermion-Z2 parafermion- as cos$.) 

The relation between the c = 2 free bosonic theory and the Potts model is not straightforward. It can 
be formulated in terms of cosets beginning with SU{3)i x SU{3)i. This conformal field theory arises 
from bosonization of critical SU{3) "spin" chains (the two factors of SU{3)i arising from left and right 
movers) and the associated boundary critical phenomena is closely related to that of the triangular lattice 
QBM problem. It will be discussed in a later paper .c3 It is natural to associate the SU{3)2 CFT with the 
diagonal SU{3) symmetry of this model. The remaining coset, SU{3)i x SU{3)i/SU{3)2, with c = 4/5, 
gives the 3-state Potts model. This is also equivalent to SU{2)3/U{1). Yi and Kane first discussed the 
non-trivial critical behavior of the the triangular lattice QBM problem using results from the 3-channel 
SU(2) Kondo problem, which is, in turn, related to the S'[/(2)3 CFT. The second factor in the conformal 
embedding, SU{3)2 can be further factorized into two C/(l) CFT's [the maximal abelian sub-algebra of 
SU{3)] and a c=6/5 coset, SU{3)2/U{lf which is naturally regarded as the Z| CFT. Since the SU{3)i 
CFT, with c=2, is actually equivalent to two free bosons, this conformal embedding essentially expresses 
2 free bosons in terms of the Potts model and the Zf CFT. One can certainly write the stress energy 
tensor for 2 free bosons as T = Ti + r2 where: 



^4 



T2 = - 



-lt(9^.f+j:- 



i(*3-*fc) 



3 



jT^k 



iEi^^.f-j:^ 



«(*3-*fc) 



j=i 



j^fc 



(5.11) 



such that the short distance expansion of Ti with T2 is trivial. Here d denotes d/dz where z — t + ix 
and T is imaginary time. Ti is a stress energy tensor with central charge ci = | of the coset SU{3)i x 
SU{3)i/SUi3)2 or SU{2)3/U{1), and similarly for T2 with C2 = | of the SU{3)2/U{1)^ coset. The 
parafermion ^1 turns out to be primary both for Ti and Ti + T2, but it is not the case for most other 
operators of interest. Some of the most crucial operators in the Potts model, like the field with conformal 
weight I , do not have any known representation in the two boson theory, neither as vertex operators nor 
generalized twist fields. 

It follows that, strictly speaking, we do not have a conformal embedding. 



1 theory: while c = 1 



situation is reminiscent of what happens at the Ising square point of the c 

and a decomposition of the c ~ 2 stress energy tensor analogous to (5.11) exists,EJ the periodic boson 



As mentioned before, the 

2 f^ 2' 
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theory cannot be considered as the product of two Ising theories; only the Z2 orbifold can. In the absence 
of an understanding of the orbifold, one would nevertheless observe that, if the torus partition function of 
the periodic boson cannot be expressed in terms of Ising model characters, boundary partition functions 
can: 



Znn = (Xo + X1/2) 

^DD - Xo +X1/2 

Znd ^xiMXo + xi/2)- (5.12) 

Again, the existence of a "boundary embedding" in our c — 2 case suggests that the situation is similar, 
and that there is some sort of orbifold version of the bosonic theor y for which there would be a genuine 



conformal embedding. Nevertheless, since the perturbation in (4.17) makes sense as a pure Potts operator 



it is natural to expect it to induce a flow purely in the Potts sector of the theory. Thus we come to the 



important conclusion that the interaction term in Eq. (4.17) is purely a Potts operator. In fact, it is 
precisely the same boundary operator which corresponds to a boundary magnetic field, hS^^i in the Potts 
model where a = 1,2,3 labels the Potts variable. The case va < then corresponds to h < 0, inducing 
a flow from free to fixed boundary conditions in the Potts model, while the case va > corresponds to 
h > 0, and should instead induce a flow from free to mixed. This verifies the handwaving arguments of 
Sec. IV relating QBM to the Potts model. Note however that this connection is only established for the 
special choice of lattice spacing a^ — An/S. 



We also consider the more general QBM model of Eq. (4.32) in which the remaining 3m point group 
symmetry is broken. Upon using the conformal embedding, we now expect that products of Potts with 
Zg ' (or Ising x tri-critical Ising) operators will appear in the Hamiltonian. 

In the remainder of this section we use the fusion technique to extract information about the boundary 
critical points of the QBM model. By considering fusion in the Potts sector we can study all the fixed 
points that occur in the 3m symmetric model. We expect that we have obtained all fixed points with this 
symmetry yielding the phase diagram of Fig. (0). On the other hand, by considering fusion in the other 
sectors we obtain a collection of additional fixed points whose properties we understand in less detail. In 
Sec. VI we discuss the integrability of some of the RG flows between fixed points and thereby confirm 
the corresponding ratios of g-factors. 

B. Fusion in the Potts sector 

The Potts sector in the first term of the amplitudes Xo + X3 a-nd Xo + X3 + 2x£/3 are exactly the 
fixed-fixed and free-free amplitudes in the Potts model. It suggests that the Dirichlet and the Neumann 
boundary conditions correspond to the free and fixed boundary conditions of the Potts model, respectively. 

In the (pure) Potts model, we can construct the free boundary state from the fixed boundary state by 
the fusion with the weight-1/8 primary operator O44 (which is absent in the bulk spectrum of the Potts.) 
In fact, the fusion with the same operator in the Potts sector in our co nfor mal embedding gives the 



Neumann b.c. from the Dirichlet b.c. The Potts sector in the first term of (5^) is the free- free amplitude 



of the Potts model, and is thus transformed to the fixed-fixed amplitude by double fusion: 

X^ + xf ^ xf/8 + xfs/s -* X^ + xf + 2xf/3, (5.13) 



where each -^ means the fusion with 1/8. On the other hand, the Potts sector of the second term of ( |5.3| ) 
is transformed as: 



X2/5 + X7/5 -" Xi/40 + X21/40 -^ 2xi/i5 + X2/5 + X7/5- (5-14) 



Thus, by double fusion with 1/8 in the Pott s sector, the Neumann-Neumann amplitude (5.4) is obtained 



from the Dirichlet-Dirichlet amplitude (5.3) 



In this fusion construction of the Neumann boundary state, the ratio of the ground-state degeneracy 
9n I go is given by the ratio of modular S'-matrix elements of the Potts model. Thus it automatically 
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agrees with the ratio of Potts degeneracy gfroc/fffixcd- The correspondence between Dirichlet/Neumann 
boundary conditions in the c = 2 theory and fixed/free boundary conditions in the Potts model is also 
consistent with the integrable field theory approach taken in Section VI. 

In a similar way, we can construct another boundary state for the c — 1 theory, which corresponds 
to the mixed boundary condition of the (pure) Potts model. In the Potts model, the mixed boundary 
condition is obtained from the fixed boundary condition by fusion with the 2/5-operator (e). Thus we 
attempt a fusion with the same operator in the Potts sector to construct a boundary state from the 



Dirichlet boundary state. The Potts sector of the first term in (5.3) is transformed as: 



Xo + X3 ^ X2/5 + X7/5 ^ Xo + X3 + X2/5 + X7/5' (5-15) 

which is the same as the transformation of fixed-fixed amplitude to mixed-mixed one in the Potts model. 
That of the second term is transformed as: 

X^/5 + xf/5 -^ X^ + Xf + XV5 + Xf/s ^ X^ + xf + 2x^/5 + 2x^/5 • (5-16) 

Thus, the amplitude for the nontrivial boundary state ( y-state) is given by 

ZYY{q) = (XoxiT + x(/2xj/2)(x^ + X.f + X^/5 + xf/s) + 

(XoxJ/5 + x(/2Xvio)(x^ + xf + 2x^/5 + 2xf/5), (5-17) 

where the argument q is omitted. By the modular transformation to the "closed string" channel, it is 
expressed as 

^yy{^) = ^TTf^KXo + XV2)(X^ + xI/2) + ^xi/wXyieKXo + X3 + ^Xys) + 



4^/3 

3- V5 

4%/3 



[(Xo + x(/2)(X^/5 + xf/io) + 2x(/i6xJ/8o](xf/5 + xf/5 + 2xf/i5), (5-18) 



where the omitted argument is now q. The ground-state degeneracy for this state is gy = 

a/(3 -I- V5)/(4\/3) — 2cos(7r/5)/v 2-^3. Again, by construction, the ratio of the degeneracy is equal 
to corresponding one in the pure Potts model. 

9N ^ m^ ^ -^. (5.19) 

gv .gmixcd 2cos7r/5 

It turns out that this Y-state is identical to the nontrivial fixed point found by Yi and KaneQ using 
a different mapping. This will be confirmed by the calculation of mobility, and also by a physical 
consideration on the relation to the Potts model. 

It has been known that the Potts model admits, fixed, mixed and free boundary conditions as con- 
formally invariant boundary conditions. Recentlytll, a new boundary condition was found in the Potts 
model and related to the mixed one by the duality transformation. The new boundary state is obtained 
by fusion with the operator O22 of dimension 1/40 from the fixed boundary state. Applying the fusion 
with the same operator to the Dirichlct-Dirichlet amplitude, we obtain a new boundary state, which we 
label W ^ for the present problem: 

Zww{q) = (XoX^ + x(/2X'V2)(x^ +X3 + X^/5 + 2xf/3 + 2xr/i5) 

+ (XoxJ/5 + Xv2Xvio)(x^ + Xf + 2x^/5 + 2x^/3 + 4xf/i5)- (5-20) 

This is another nontrivial boundary state, which corresponds to the Potts "new" boundary state. The 

ground-state degeneracy of the WK-state is given by gw = 2 cos (7r/5)Y \/3/2. 

Thus, there is a corresponding boundary state in the present problem for each boundary state in the 
Potts model. The ratios of g-factors between the boundary states are identical to the corresponding ratios 
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in the Potts model, by construction. Moreover, more boundary states can be constructed in the Potts 
model by forming superposition of several boundary states. Such boundary states contains dimension-zero 
boundary operator(s) other than identity. Although they are often unphysical due to their instability, 
they can be relevant for some physical situations with a first-order transition. 

In the present problem, the potential minima form a triangular lattice isomorphic to T for va > 0. 
However, for negative va, the potential minima form a hexagonal lattice, which has two points within 
each unit cell of T. In the limit va — > — oo, the corresponding boundary condition is Dirichlet, but 
the boundary field is pinned to one of two inequivalent minima. This can be represented by "Double 
Dirichlet" boundary state, which is the superposition of two Dirichlet boundary states for the two minima. 
Considering the correspondence Dirichlet ^^ fixed and Neumann ^^ free, it would be natural if the double 
Dirichlet corresponds to the sum of two fixed, say A and B boundary states in the Potts model. Indeed, 
the double Dirichlet amplitude can be expressed as follows: 

ZoD^DDiq) = {xixl + x[,2xl/2){'^Xo + 2xr + Sxf/g) + {xixl,^ + Xv2Xvio)(2xf/i5 + Sxf/s), (5-21) 

where the Potts part of the first term is the amplitude for the superposed state A + B. 

To summarize, based on the conformal embedding, we find several boundary states for the c = 2 
problem constructed by fusion in Potts sector as follows. There is a c = 2 boundary state corresponding 
to each Potts boundary state. 

c = 2 model Potts model ^-factor (for c = 2) 



Dirichlet Fixed {A) l/VVl 



Neumann Free a/\/3/2 

Y(i-Kane) Mixed (AS) 2cos(7r/5)/V273 



Double Dirichlet A + B \/ 2/\/3 

W New 2cos(7r/5)y^V3/2 

C. Fusion in Ising/Tricritical Ising sector 

It is also possible to construct other boundary states by fusion in Ising or tricritical-Ising sectors. 
The U-state is obtained from D-state by fusion with cither 1/16 operator of Ising or 7/16 operator of 
Tricritical Ising. Its diagonal partition function is given by 

Zuu{q) = (Xo + X(/2)[(X^ + xl/2)i.X^ + xf ) + (Xvio + xl/5)X2/5l (5-22) 

with the degeneracy gu — 1/3^'*. Similarly, S-state is obtained by fusion with the same operators from 
the N-state. 

Zssiq) = (X^ + X(/2)[(X^ + Xl/2){X^ + X^ + 2xr/3) + (xf/io + xl,^){'ix'^/i, + xf/5)], (5.23) 

with the degeneracy gs — 3^'*. It turns out that these state represents mixtures of Dirichlet and 
Neumann: imposing Dirichlet in one component and Neumann in the other. 

On the other hand, another boundary state, the T-state is obtained from D-state by fusion with 3/80 
operator of Tricritical Ising. The diagonal amplitude is 

ZMq) = (Xo + X(/2)[(X^ + Xf/io + xJ/5 + Xl/2)XI + (x^ + 2x^/10 + 2x^/5 + xl/2)x.l (5.24) 

with the degeneracy gT = (7 + 3\/5/6)^''''. Yet another state, R, is found by fusion with the same operator 
from N-state. 
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ZB.R{q) = (xo + x(/2)[(x^ + xvio + xJ/5 + xj/2)(x/ + xv + x^-t) 

Hxl + 2xvio + 2x^/5 + Xj/2)(X. + X.t + X.)], (5.25) 

which has the degeneracy qr — [3(7 + 3-\/5)/2]^/^. The calculation of the mobility, which will be given 
later, suggests that they are mixtures of Y and W, namely taking the dual of one component of the 
bosons in Y or W boundary states. 

D. Fusion in the Z^ sector 

(5) 

We can apply similar lines of arguments as in the previous section to the Zg sector in the conformal 
embedding of the c = 2 model. As we will see, we obtain new boundary states as well as those found 
using the Potts sector. Starting from the Dirichlct boundary state, we have tried fusion with all primary 

(5) 

fields in the Z^ model. As a result, we found new boundary states, besides the same ones obtained 
by Potts fusion. For simplicity, here we focus on the boundary states with single dimension-0 boundary 
operator. 

A new boundary state F is obtained by fusion with Z^^' operators of dimension 1/9, 7/9 or 13/9 from 
the D boundary state: 

ZpF = Xi{xl + '2xl + ixl,2) 

+X.(2X3/5+X^/5 + 3x?/io)- (5.26) 

This turns out to represent free diffusion under a magnetic field,E3. While a detailed discussion of this 
identification will be given in a later publication, we will see a signature of the magnetic field, namely 
the Hall effect, in the next subsection. 

Another one, X is obtained from the D by fusion with an operator of dimension 2/45, 17/45 or 32/45. 
The amplitude in the open string channel is given by 

Zxx{q) = Xi{xl + 2x3/5 + X8/5 + 2X2 + 3x^/2 + 3xvio) 

+X.(X^ + 4x^/5 + 2x^/5 + 2x^ + 3x?/2 + 6x?/io) (5-27) 

By construction, the ground-state degeneracy gx satisfies 

9D 5° 

(5) 

where S is the matrix modular transformations of the Zg model and and c represents the identity 

operator and the operator used in the fusion. S^ — y (5 — a/5)/10/6 and S*" — y (5 + a/5)/10/3 gives 
gx = 4cos (7r/5)5£) = 2gy. Fusion with an operator of dimension 3/5 or 8/5 gives the same F-state 
previously obtained by fusion in Potts sector. 

On the other hand, starting from the Neumann boundary state, we obtain two additional boundary 
states V and Z by fusion in the Z^ sector. No other states are found by fusion in Z^ on other known 

(5) 

states. To summarize, there are four new boundary states i^,X, V and Z obtained by fusion in Z^ 
sector. As we will see later, these four states exhibits finite Hall mobility, and are presumably related to 
RG fixed points of QBM under a magnetic field. 

E. Mobility and nature of the boundary states 

Now let us calculate the mobility in the QBM problem. It helps us to identify the physical natures of 
the constructed boundary states. The mobility is defined as the linear response coefficient of the velocity 
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of the particle to the external force. The external F acting on the particle is represented by the additional 
term to the Hamiltonian, —F ■ Q. The linear response of the velocity of the particle to this force is given 
by the Kubo formula 

-*">- f. /,,(*32)<3.„)). (5.28) 



^ dt ' J ' dt 

Thus the static mobility tensor can be expressedlH as 

,..,= lim^(^H§^(-.)). (5.29) 

i^^o u dt dt 

Using the weak corrugation formulation we can identify Q with 0(0), and the velocity operator dQ/dt 

with d(f>{0)/dt, which corresponds to the current operators. Thus the mobility is deduced from the 
two-point boundary correlation function of the current ja = Ja + Ja (o- — x,y), where Ja and Ja 
are a-component of the holomorphic and antiholomorphic currents,, pThe mobility corresponds to the 
conductance in the problem of a tunneling between quantum wireua, which can be also mapped to 
essentially the same boundary CFT problem. Note that we now have a two-component boson, and thus 
the current has two components. 

The time dependence of the correlation function is basically insensitive to the boundary, although 
some coefficients do depend on the boundary state. As in Ref. ( ^2|), the JaJb and JaJb correlations are 
insensitive to the boundary conditions, and can be fixed by an appropriate normalization of the current 
operator. On the other hand, while the Ja ■ Jb correlation obeys always the same power law, thepamplitude 
depends on the boundary condition. The amplitude Aab is determined by the boundary statec3 as 

_ {Q\Ja,Jbs\B) 
Aab- ^Q|^^ , (5.30) 

where |0) is the ground state, \B) is the boundary state under consideration, and Ja,i {Jb.i) is the level- 
1 annihilation part of the holomorphic (antiholomorphic) current. By a suitable normalization of the 
current operator, the mobility /iat can be written as 

Oab + Aab /c oi ^ 

^^ab = ^ . (5.31) 

Since the N boundary condition corresponds to the delocalized phase, we normalize the mobility in the 
N phase to be unity, fixing A^ = 5ab for the Neumann boundary condition. On the other hand, the D 
boundary condition corresponds to the localized phase and the mobility should vanish, thus A^^ = —Sab- 
In our conformal embedding approach, the two level- 1 states are identified as 

J.,-i|0) = |0)/|3/5)t|2/5)p, (5.32) 

J^,_i|0) = |1/2)/|1/10)t|2/5)p, (5.33) 

since they are the only two combination of Virasoro primaries with weight 1 . 

The coefficients Aab of these states in the boundary state are changed by fusion. The changes are 
given by the modular ^-matrix elements. Thus we can calculate the coefficients and consequently the 
mobility of the new boundary state constructed by fusion. Since the both states contains the same Potts 
primary state \2/5)p, fusion in the Potts sector gives Aab — A6ab- Thus the mobility is also diagonal and 
isotropic: fiab = fJ-Sab- 

The TV-state is obtained by fusion with the operator of dimension 1/8 in the Potts sector from the 
D-state, and the amplitudes is given by 



2/5 



A^ Stl 50 



AD 2/5 ^0 



= -1, (5.34) 
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where S is the matrix of modular transformations for the Potts modeL This is indeed consistent with 
the physical consideration fi^ — and /i^ — 1. 

Similarly, the Y-state is obtained by fusion with 2/5-operator in the Potts sector and thus the amplitude 
for the Y-state Ay is given by 

The mobility /i is then 

A*^ = (5-\/5)/4 = 2sin2 7r/5. (5.36) 

This agrees with the result obtained for the nontrivial fixed point by Yi and Kane,Ll implying that our 
y-state is identical to their fixed point obtained by a different mapping. 
The W-state is obtained by fusion with 1/40-operator in the Potts sector: 

A^ ^Sl/% so _ 3-V5 

Thus /i^ = (VS — l)/4. This is equal to 1 — /i^; this is related to the duality between the F-state and 
VF-state. 

Next let us consider the fusion in the fsing/tricritical fsing sector. Since the two components of the 
currents are related to the Ising/tricritical Ising sector in an asymmetric manner, in general the amplitude 
Aab is asymmetric {A^x 7^ Ayy) although it is diagonal {A^y = 0). For example, the U-state is obtained 
from D-state by fusion with the 1/16 primary in the Ising sector. Thus, the coefficients are given by 

Axx — co "co — ^xx = ^1' (5.38) 

"^0 "^1/16 
1/2 
aU _ 1/16 ■S'o .£) _ -. ._ oq^ 

VV ~ 0I/2 QO VV ~ ' \o.ojj 

So -^i/ie 

where S represents the matrix of modular transformations of the Ising model. As a result, the mobility 
in the U-state is anisotropic: /i^^ = 0, fiyy = 1 {^^y = 0). This represents perfect mobility in j/-direction 
and complete localization in the orthogonal a;-direction. This would naturally correspond to a "mixed" 
b.c. of D and N. Indeed, this is the case as we demonstrate in the following. 

A similar calculation gives the mobility for the S-state as /xf^ = 1, ^y = (/if„ = 0). One might think 
that U and S are equivalent upon a space rotation, because their mobihties are the same if we exchange 
X and y directions. However, they are not equivalent because the underlying lattice F is not invariant 
under tt/2 rotations. The inequivalence can be also seen in the following construction of the mixed D/N 
boundary state. 



As discussed in Section IIIB, the mixed D/N state can be constructed from (2.63) with an orthogonal 
matrix R. Let us choose 

This would correspond to D (localized) in x-direction and N (free diffusion) in y, as in the U-state. Then 
the allowed zero modes are given by the integer-coefficient linear combinations of 



[Ui,Vi\ 



(-^,0),(0,0) 



(5.41) 



[^2,1/2] - [(0,0),(0,2V^)]. (5.42) 
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Thus, the new lattice f introduced in Section IIIB is a rectangular lattice with the unit cell of the size 
— i= X -y^. The diagonal amphtudc for this state is given by 



2 



^('^)-(^) El''^^'^'- (5.43) 

This actually agrees with the diagonal amplitude Zuu{q), implying that the U-state is identical to the 

constructed mixed D/N state. Of course, the ^-factor also agrees: g = J TrVaiT) = 3^^/"* = gu- 
On the other hand, the S-state would correspond to 

R-{1\)- (5.44) 



In this case, the bases for the allowed zero-modes is given by 

[uuvi] = (0,^),(0,0) , (5.45) 

[U2,i?2] = f(0,0),(2\/3^,0)l . (5.46) 



This is inequivalent to the previous case ( ^.40| ), because the lattice T is not invariant under a rotation by 
7r/2. Now the new lattice T is again a rectangular lattice, but with a different size I/a/tt x -^S/tt. This 



gives g = •W7rVo(f ) = S^''* — gs- Moreover, its diagonal partition function agrees with Zss, implying 
that the S-state is in fact the mixed D/N state. 

We can construct other mixed D/N states, by considering other (infinitely many) possible rotation 
matrices R. The fusion construction based on the present conformal embedding does not give those 
states. It does not mean that these states do not exist. It rather means that, only a part of the possible 
b.c.'s can be accessible by the conformal embedding which reduces the number of conformal towers 
effectively. In general, we do not know how far we can reach by a particular conformal embedding. From 
physical considerations, however, we think that the fixed points appearing in the QBM model introduced 
in Section IV have been exhausted. 

The T-state is obtained from the D-state by fusion with 3/80 primary in the tricritical Ising sector. 
The coefficients are given by 

■ T ^ ■^3/80 5'g jj ^ 3- %/5 
^^ o3/5 C-O ^xx r, ' 

Sq' '^3/80 ^ 

■ T ^ '^3/80 5'g jj ^ -3 + \/5 

yy oi/io qo yy o 

Sq' -^3/80 ^ 

Thus the mobility is given by /i!^^ — fi^ and /i?^ ~ /i^ = 1 — /jX . The 
Y and W boundary states. 

By a similar calculation, we find the mobility in the R-state as /z^^ 
mobility in the R-state is equivalent to the T-state after tt/2 rotation in the ccy-plane. They are again 
not equivalent, presumably reflecting the fact that the lattice T is not invariant. T- and R- states are 
apparently given by "taking dual" of Y or W for only one component of the boson. However, we do not 
know how they can be realized in text of QBM. Possibly they correspond to some anisotropic QBM. Below 
we summarize the properties of all the boundary states obtained by fusion in Potts, Ising or tricritical 
Ising sectors. 









(5.48) 


tatc 


seems 


to be 


a mixture of 


M^ 


f^yy - 


-^^y. 


Namely, the 



48 



Boundary state g-factor fx^x fJ'yy 



D l/\/2V3 



N V\/3/2 1 1 

Y 2cos(7r/5)/ V2V3 (5 - \/5)/4 (5 - V5)/4 

W 2cos(7r/5)y^V3/2 (-l + V5)/4 (-1 + V5)/4 

R [3(7 + 3\/5)/2]i/4 (-l + V5)/4 (5 - V5)/4 

s 31/4 1 

T (7 + 3V5/6)i/4 (5-V5)/4 (-1 + V5)/4 

U I/31/'* 1 

In the embedding using Z3 , the level-1 "current" states corresponds to the product state ||)5||)p, 

where ||)5 is the primary state of the Zg theory with weight 3/5. In fact, there are two such primary 

states in the Z^ ' theory, which arc complex conjugates. Let us denote one of them as || ). In this 
approach, the natural identification would be 

i(J.,-i+^J,,-i)|0)-|)5|^>p, (5.49) 

i(J,,_i-*J,,_i)|0) = |^*)5|^)p. (5.50) 

Assuming this correspondence, for a boundary state constructed from the D-state by fusion with the Z3 
operator 0, we find 

Aab = -SabRcA - eafelmA, (5.51) 

where eab is the antisymmetric tensor e^y = —£yx = 1, ^xx = ^yy = 0. Here A is defined as 

c.3/5 cO 

O ± On, 

A^4rj^, (5.52) 

q3/5 qo ^ ' 

where 5^ is the matrix of modular transformations of the Z^ theory. Thus, when the matrix of mod- 
ular transformations acquires an imaginary part, there is a non-vanishing off-diagonal mobility. This 
corresponds to the Hall effect. Physically, such an effect is expected in QBM under a magnetic field. 
The F-state is obtained by fusion with the 8/5-operator in the Z3 ' sector. Thus 

._^_^_ 3-a/5 
^'STSI,,-- 2 • ^'-'^^ 

This gives the same result as we have obtained by fusion in the Potts sector. The off-diagonal mobility 
vanishes in this case. 

On the other hand, a similar calculation on the F-state reads 

. _ -^i/g 5*0 -l-\/3^ 
' " V^'W, = 2 ■ (^-^^^ 

Oq 1/9 

This complex amplitude gives a non- vanishing off-diagonal mobility: 

t^-xx — ^J-yy — .J fJ'xy — fJ-yx . • (_0.00j 
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By a fusion with the conjugate operators, we can also obtain the boundary state with opposite sign of 
the off-diagonal mobility. 

The other three states X,V and Z also exhibit the non-vanishing off-diagonal mobility. Thus they 
should also correspond to a critical behavior under a magnetic field. The mobility in these states is 
summarized as follows. 



Bo 
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state 


g-factor 
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f^xy 
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_|_y3 
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[f -f2V5]i/4 


1 + ^5 
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_i_3V3-VT5 
-^ 8 
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F. Structure of the non-trivial boundary state 



We have obtained the explicit expression of the amplitude for the boundary state Y. ft gives some 
insight into the structure of the nontrivial boundary state. 

In general, a conformally invariant boundary state should be a linear combination of Ishibashi states. 
Each Ishibashi state is constructed f rom a spinless primary field that appears in the bulk theory. Thus, 
the closed string channel amplitude (5.18) must be consistent with the bulk spectrum. 

In c = 1 free boson theory with a generic compactification radius, each bosonic Fock space build on 
a vacuum corresponds to an irreducible representation of the Virasoro algebra (except for zero winding- 
number sector.) However, in c = 2, a bosonic Fock space is always reducible and contains an infinite 
number of irreducible representations of the Virasoro algebra. Thus, there are infinitely more Ishibashi 
states other than those built on a bosonic vacuum. Nevertheless, the c = 2 Dirichlet and Neumann 
boundary states have simple expressions in terms of bosons, and their (self-) amplitudes are written as 
Sft, Q^ I'ni^iY ■ Roughly speaking, the factor l/rj{qY means the boundary state is made of whole bosonic 
Fock space. 

However, the non-trivial boundary state is not made of such "bosonic" boundary states. This can 
be seen as follows. Expressing the diagonal amplitude (5.18) as a sum of the character if' /r]{q)'^ of the 
Heisenberg algebra. 



ZvYiq) = 



9Y 

v{q) 



1 + (21 - 9^5)9^/'^ + ^(1-^5^1/2 + (16 - 6V5)g2/3 



-f(5 - 3V5)q + 9(3 - V5)g^/« + 3(7 - 3V5)f^^ + (5 - 3^5) f/"" + 



(5.56) 



This involves negative coefficients at least for q and (f''^. While any given partition function can be 
written as an infinite sum of the Heisenberg character, the coefficients reveal the nature of the boundary 
state. If the boundary state is a linear combination of the bosonic boundary states, all the coefficients 
must be positive. 

Thus , this is a highly nontrivial boundary state which cannot be constructed by a generalization such 
as eq. (2.63) or as in Ref. ( 30) of the Dirichlet or Neumann. To our knowledge, this is the first proof 
that such a non-bosonic boundary state does exist in a free boson field theory. Similar non-triviality can 
be shown for W,R,T,X and V states. 

On the other hand, since the boundary state is a linear combination of Ishibashi states, the diagonal 
amplitude must be a linear combination of c = 2 Virasoro characters with positive coefficients. In fact, 
the amplitude ( 5.18| ) is expressed as 



4V3Zyy(q) = (3- 



V5)xl 



6(3 - V5) 



Xl/6 



18x?/2 



(18-2y5)x^/3 



(5.57) 
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with positive coefHcients, where the c = 2 Virasoro characters for weight h are given by Xh — 
^/»-i/i2/-Qoo^^(l _ ^«) ioT h > and xo = (q'^^^^ - ?"/'')/ 11^=1(1 " <?") ^r h = 0. Moreover, 
all the primary weights correspond to those of spinless fields in the bulk spectrum. These are guaranteed 
by the fusion construction. 

Thus this nontrivial boundary state is certainly a consistent boundary state for the c = 2 free boson 
theory, although it appears to be very complicated in terms of the bosons. A possible direction to extend 
the construction of the nontrivial state for different lattices F* is to write down the boundary state in 
the bosonic language and then guess a generalization of the boundary state for other lattices. However, 
so far we have been unable to develop this approach. 



VI. INTEGRABLE FIELD THEORY ANALYSIS 
A. Integrable flows away from the perfect mobility (Neumann) fixed point 



Here we consider the Hamiltonian of Eq. (4.32|) for small vi, V2, v^. This model turns out to be 



integrable for general values of vi, V2, v^ as discussed in detail in Ref. ( ^). This integrability is, in a 
distant sense, related to the existence of the conformal embedding; more precisely, it is the underlying 



parafermionic algebra satisfied by the verteXj-Qperators of (4.17) at that particular radius that ensures 
the existence of non local conserved currents.a 



For general values of Vi, V2, v^ the scattering theory describing (4.32) is simply made up of three left 
moving and three right moving massless particles with mass parameters mi, rM2, 1^3 , energy and momenta 
being parametrized by rapidities e = ±p = nije^ . Scattering between left and right movers is trivial since 
the theory is scale invariant in the bulk, while scattering among left or among right movers is described by 
pure CDD factors: Sk,k±i = itanh (| — ■^), ^i.a = 1. Only the first particle then scatters non trivially 
on the impurity, with a boundary R matrix again given by a CDD factor, R — itanh (^^^^ — ■^). Here, 
6b parametrizes the impurity energy scale Tg oc e^^ ex vf . The exact dependence of Tb, as well as the 
dependence of the mass parameters upon the v^'s are complicated and will not be needed in the following, 
except for special cases. 

The main result of the integrable approach is that the impurity free energy can be computed via the 
thermodynamic Bethe ansatz (TBA).E3 If we parametrize the filling fractions of the particles by pseudo 
energies ej as ij = t-jt 1 the latter are solutions of a set of integral equations 

" = ""'' - ^S ».. / f ;s^ '» (1 + =-*"'"^) . (") 

Here, j, k take values 1, 2, 3 and can be represented as the nodes of an ^3 Dynkin diagram with incidence 
matrix Njk (crossed nodes indicate the existence of a source term in the TEA equations) 

mi ?7i2 TO3 




The impurity free energy reads then 



Formula (x2) gives us quick access to the change of boundary entropies between the UV fixed point 
(Neumann boundary conditions, T /Tb = cxd) and the possible infrared (IR) fixed points [T/Tb ~ 0) our 
model can fiow to. Rather than entropies, we will use the ground state degeneracy, g = e'^/"^. The change 
of g-factor is then neatly expressed as 
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guv _(l + e ^ ^M . (6.3) 



gm V 1 + e-<^i(°°)/'^ 

When 9 -^ — cxd, the source terms just drop from the equations, the e's go to constants that obey the 
system, setting Xj — e~'j(~°°)/'^, 



cj =n'^^+^'=) 



Wjfc/2 



This is solved right away by xi = x^ = 2, 2:2 = 3. 

Setting similarly yj = e~'^j(°°)/^, the system obeyed by the j/j's depends on the source terms: different 
possibilities can arise. If toi ^ 0, the source term for ei diverges, and therefore j/i = 0. If mi — and 
m,2 7^ 0, in the IR the system is simply ei = so j/i = 1. Finally, if 7711 = and m2 = 0, the IR system is 

Z/i = y2 = (1 + yi)^^^ 

with solution j/i = ^2^ ■ ^'^ thns obtain the possible ratios of degeneracies 

giR 

9UV U a -in 
= \/ o , mi = 0, m2 7^ 

giR. V 

= I^sfx, mi = 7712 = 6.4 



with 2 cos r. — M 2 



3+75 



We obtained the first ratio in Sees. Ill and IV when discussing the flow from N boundary conditions to 
D boundary con dition s (i.e. freely diffusing to localized) in the case vi = V2 = v^ > 0. For general Vi, the 
potential of Eq. ( 4.32| ) has a unique minimum (within each unit cell) and the TBA equations in the first 



case describe the flow from N to D in the general case where th e th ree-fold symmetry is broken. This is 



confirmed by the analysis of the central charge associated with ( |6.l| ) which gives generically c = 2, hence 
indicating a flow within the whole two boson theory, (i.e. not purely within the Potts sector.) The third 
ratio was also obtained previously when we discussed the flow from N to Y in the case vi — V2 — v^ < 0. 
Since this ratio is obtained only when mi = m2 = 0, it follows that the flow from N to Y occurs only in 
the three-fold symmetric case. Meanwhile, the analysis of the central charge associated with (g]j) gives 
only c = I , indicating that the flow takes place within the Potts sector of the theory only. 

Finally the second ratio is somewhat trivial: it corresponds to a flow from Neumann to Dirichlet 
boundary conditions for the field 4>i while the field 02 remains at the Neumann fixed point. This 
corresponds to the U boundary state discussed in Sec. IV. This is obtained by s etting V2 — vj, — 0, while 



vi can take either sign. The analysis of the central charge associated with (6T) gives c = 1, confirming 
that the flow takes place in the c = 1 sector only. 

B. W3 integrable flow^s 

Integrable flows can be used to explore more thoroughly the phase space and the possible fixed points 
in our problem. 

The integrability of the flow from free to fixed or mixed boundary conditions in the Potts model, and 
therefore, from Neumann to Dirichlet or the Y fixed point in the c = 2 theory, is ultimately related 
with the fact that the Z^ parafermionic theory (the three state Potts model) can be represented as the 
coset -Tjff) ■ The field ^1 is then the "adjoint" operator in this construction.]^ We observe that this 



^In a general coset construction, one usually calls adjointH the operator obtained by taking the identity repre- 
sentation for the algebras in the numerator, and the adjoint for the algebra in the denominator. 
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theory can also be represented as the coset — sij(3) ■ 



The adjoint operator for this coset now has 



dimension A = ^, and physically corresponds to the energy operator of the Potts model. It can be shown 
to define an integrable perturbation that preserves the W3 symmetry: more precisely, there is a natural 
deformation of the current W3 that is conserved, at least perturbatively, away from the UV fixed point, 
if this fixed point itself preserves W3 symmetry. It is natural to assume that this extends all the way to 
the IR fixed point, which therefore should have W3 symmetry too. 

In our problem, the operator with A = | is present in the spectrum for mixed boundary conditions 
in the Potts model, or, equivalently, at the Y fixed point in the c = 2 theory. In the Potts model it 
corresponds to applying a magnetic field at the mixed (AB) fixed point which breaks the remaining Z2 
symmetry and favors the A state. This corresponds to a flow from mixed to fixed. In the QBM problem 
we expect this operator to produce an RG fiow from Y to localized fixed points; if Y corresponds to the 
particle being on A and B sub-lattice then the fiow is to a state where the particle is localized on an A 
site. This RG fiow can also be studied using integrability. 



For the more general coset 



SU{3)ixSU{3)k 
S(7(3)fc4 



perturbed by the operator of weight A = 1 — j^ , which 



e ope 
.n@, 



preserves W3 symmetry, S matrices and R matrices are knownErl, and the impurity free energy can be 
computed once again using the thermodynamic Bethe ansatz. It reads 



F 



-T 



de 

2^ 



Gi( 



,)ln 1 



'Cii/T 



G, 



lln 1 



^e2i/T 



(6.5) 



Here, as before, Tb is an impurity energy scale, Tg oc e^^ , Gi and G2 are known kernels, and the e are 
pseudo energies, solutions of the equations: 



ey = <5jime^ -TJ2 ^^J■^kGl * In (l + e"'"=^^) -^Y, ^^.0^2 * hi (l 



.<^ij/T 



(6.6) 



Here, ij are li ne a nd column labels on the following diagram (the cross on the (1, 1) node stands for the 
mass term in (k.d 



o — o- 



O 



-o- 



k 

-o 



o — o 



i — 1,2, j — I, . . . ,k, and Nij^ki its incidence matrix. 

We will not need the exact form of the kernels in what follows, and shall concentrate on the ground 
state degeneracies in the UV and in the IR. From the previous formulas, we find the simple result (the 
Hji all vanish here) 



^^ = i[ln(l + xn) + ln(l + X2i) 
9iR 2 



(6.7) 



where xa = e "^'^'^ is solution of the 



-00 limit of the TEA system: 



=n(i+--f--^^n(i+;^ 



(6.8) 



For the case fc = 1 of interest here, the system reduces to x = xn — X21, x ~ /^ or a; = 2cos^. 
Hence 



9uv r, Ti" 
= 2cos— . 

9iR 5 



(6.9) 
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This is the ratio ffniixed/^fixed ^^ ^^^ Potts model, or the ratio gy/go for the two boson system. We 
thus expect that perturbation by the operator A = | does not lead to a new fixed point, but induces 
a flow from Y to D in the two boson system. Let us stress that this is a rather abstract statement: 
although we know the operator with A = | is present in the spectrum (because of the partition functions 
analysis), we are not able to represent it in terms of bosonic operators, since we do not have an entirely 
clear picture of what the Y boundary conditions mean for the bosons. 



VII. GENERALIZATION TO HIGHER DIMENSIONS AND RELATIONS WITH THE 

KONDO MODEL 



A. Generalization to higher dimensions 



Many of our results can be generalized to the case of n — 1 bosons, corresponding to QBM on an 
(n — l)-dimensional lattice. Yi and Kane found a generalization of the Y fixed point for all n > 3 
using a mapping onto the n-channel SU{2), overscreened Kondo model. They also observed that the 
3-dimensional {n — 4) case corresponds to QBM on a diamond lattice. The conformal embedding is now 
naturally understood in terms of an SU{n) "spin" chain. We may decompose the product of left and 
right moving SU{n)i factors into SU{n)2 x a Z„ parafermion. We may then divide out the maximal 
abelian sub-algebra, C/(l)"~^ from SU{n)2. Taking into account that SU{n)i is equivalent to n — 1 free 
bosons, we see that we have constructed a conformal embedding: 



n — 1 = 



2(n-l) n(n-l) 



(7.1) 



The first component of the embedding is the theory of Zn parafermions, that can also be considered as 



the coset 



SU{2)„ 



with 



_ 2(«-l) 



The fundamental parafermion in this theory has weight - — -, and can 



U{1) ' "'"" " n+2 

be bosonized using a system of n bosons satisfying (for the chiral components) 

($,(z)$,(w;)) = -2^lniz-w) 
($,(z)$j(w;)) = |ln(z-u;), i 7^ j, 



(7.2) 



as 



*i = E' 



,i*fc 



(7.3) 



fc=i 



The embedding then follows from the decomposition T = Ti +T2, where 



Ti 



T2 = 



'.tm.f+E' 



,i(*j-*fc) 



i=i 



j^k 



^(a<i>,)^-^e^(*^-*^) 



J=l 



j¥fc 



(7.4) 



such that the short distance expansion of Ti with T2 is trivial, Ti is a stress energy tensor with central 
charge ci = 2^ of the coset SU{n)i x SU{n)i/SU{n)2 or 5t/(2)„/t/(l), and similarly for T2 with 

C2 — n+2 of the SU{n)2/U{l)'^^^ coset. T — T1+T2 of course is the stress energy tensor of the initial 
(n — 1) boson theory. 

It turns out that the n — 1 boson theory with boundary perturbation 



^fjCOS$j 



(7.5) 
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is integrablea, the particular case where all the Vj 's are equal corresponding again to a perturbation of 
the form '^i + \1/|. The impurity free energy reads, generalizing (6.2) 

d0 1 



^™P^ ^ ' 2^cosh(0-0B: 



■ln(l + e-'="W/^), (7.6) 



where the e satisfy the same TBA equations as ( p.l\ ), but the incidence matrix is the one of the following 
diagram 

n — 1 
1 2 n-3 "^ 

— (8) 0- - -O — (8) ^ - 2 



Here again, Tb oc w", and the tuj are dimensionless numbers functions of the couplings Vj 
In the UV, the system of equations for the x's is solved easily with Xj = (j + 1)^ — 1, j 
Xn-i — Xn — n — 1. If ?TT,„ docs not vanish, one has yn — 0, and then we obtain the ratio 

9uv 1- 
= -sjn. 

9iR 



(7.7) 



This generalizes the first ratio in (6.4), and corresponds to a flow from Neumann to Dirichlet boundary 
conditions in the n — 1 bosons problem, where at the D fixed point, the field lies on a ri — 1 dimensional 
lattice generalizing F* . 

When 77i„ vanishes, the ratio of degeneracies can take various values depending on which of the other 
TOj's vanish. A particularly interesting case is when all nij but rrin-i vanish. This corresponds to a 
perturbation with all aj ~ 0, j =^ 1, and ai > 0. In that case, one gets a ratio 

guv -sfn 



giR 2 cos ; 



(7.5 



generalizing the case n ~ '^ we denote the corresponding fixed point by Y^- These fixed points were 
discovered by Yi and Kane. El 

Let us discuss the other cases. If to„_i still does not vanish, and in addition some of the other tti^'s 

don't vanish (j < n — 2), the ratio is of the form ^^iiii — _ — ^"^ indicating a flow to a fixed point 

9 1 li ^ cos ^_|_2 i^ 

encountered previously for a lower value of n, Yn-k- When m„_i does vanish, while some of the other 
mj's don't {j <n— 2), one gets a ratio of the form Siiv_ — lji_^ xhis corresponds to an IR fixed point 

where k components of the field have D boundary conditions, the other ones having N. We see therefore 
that, as in the n = 3 case, no other fixed point is reached, besides the various D, N combinations, and 
the Yk fixed points. 

As before, the perturbation of N boundary conditions by the adjoint operator in the coset ^Ls" does 
not preserve the Wn symmetry. A perturbation preserving this symmetry is the one with the adjoint 
in the coset picture — siUn) ' ^^^ perturbing field now has dimension A = 1 — ■:^^, and the 



perturbation is integrable. The boundary free energy reads as in (6.5), with now a sum over n—1 pseudo 



energies. The TBA diagram is as in the SU{3) case, but the "base" is the Dynkin diagram of SU{n), ie 

T{n)ixSl 
SU(n)k-i 



An-1 instead of A2. For more general cosets — suT) — ^^i the diagram looks as follows 
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o — o- 



■ ■— o- 



k 

o 



n — 1 



o — o- 



-o — o 



® — o- 



...— o — O 1 



with the obvious generalization for the equations (6.6) and (6.5) that the row labels now run from 1 to 
n — \. In particular the ratio of boundary entropies in the UV and IR now reads 



guv 

giR 



-. n— 1 

1=1 



ln(l + Xii) . 



For k ~ \ we have 



Xil 



n^ 



Xjl 



-W.i,ji/2 



(7.9) 



(7.10) 



The solution of this system is 



/ sinj7r/n+2 \ '' 
\ sin7r/n+2 J 

guv _ gY_ 
giR go 



— 2 cos 



1. By some simple manipulations one finds then 

7r 

^- (7.11) 



This corresponds to a flow from y„ to D in the n— 1 bosons language, generalizing ( |6.9[ ). 

As in the n — S case, the key ingredient to understand the phase diagram is the set of boundary 
conditions and flows in the Z„ parafermioBS. theory SU{2)n/U{l). These theories have known lattice 
model or quantum spin chains realizationsEj, and it is possible to show, using results in Ref. ( |49| ), 
that the N and D fixed points correspond to free and fixed boundary conditions, while the operator 
^1 + ^1 is the most relevant order parameter, and corresponds again to a longitudinal boundary field. 
The generalization of the mixed boundary conditions for the {Z3) Potts model to the Z„ case has not 
been investigated, to our knowledge. 



B. Relation to the Kondo problem 

The emergence of parafermions is also closely related with the n-channel Kondo problem. Recalltj 
that this problem involves originally 2n Dirac fermions, which can be bosonized in terms of the current 
algebras SU{2)n x SU(n)2 x U{1). Only the spin currents interact with the boundary spin, leading to a 
fiow entirely within the SU{2)n sector, from the weak coupling fixed point in the UV to the non Fermi 
liquid Kondo fixed point in the IR. The ratio of degeneracy factors is guv/giR — 2/(2 cos ;^)- 

As is well known, the strong coupling Kondo fixed point does not depend on the anisotropy, and can 
be reached as well starting from an anisotropic Kondo interaction. There is a particularly simple choice 
of this anisotropy which corresponds to the Toulouse limit in the n = 1 case, or the Emery Kivelson 
limit in the n = 2 case, where an additional C/(l) decouples, and the flow takes place entirely within the 
SU{2)n/U{l) sector, ie the parafermionic theory. The perturbation Hamiltonian then reads 



Hi, 



V (*lcr" 



*icr^ 



(7.12) 



One can then argue from the integrable analysis that the IR fixed point - the strong coupling Kondo 
fixed point - coincides with the Yn fixed point identified previously, the generalization of the "mixed" 
boundary conditions of the Potts model to the Zn case. The UV fixed point - the weakly coupled Kondo 
fixed point - does not have a simple interpretation in terms of the Z„ variables. For the case n = 3, it 
coincides with the "localized on B or C" fixed point discussed in section IV. This correspondance was 
used extensively by Yi and Kane 13 
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C. Another relation to the Kendo problem: a remark on the marginal case 

Returning to the two boson model {n = 3), we note that the case a^ = ^ can also be studied by 
integrability techniques, using an unexpected mapping onto the four-channel Kondo model. To do this, 
recall that it is possible to represent the currents of the SU{2)4 algebra (that has c = 2) with two bosons 
as follows ;n 



Jx = (^x cos(V27r(/)iL + \/Q'k(J)2l) 
Jy = ay cos(%/27r^iL - \/67r02L) 
J,= a,cos{VS^^iL). (7.13) 

A similar representation was used by Fabrizio and GogolincZI recently; these authors however did not 
explicitly consider the a operators in their representation (cocycles), which will play a crucial role in 
what follows, though it doesn't affect their results. The four channel anisotropic Kondo problem has a 
boundary interaction 

H,nt = J± [txJx{0) + Tyjy{0)] + J//r, J,(0), (7.14) 

and is integrable for any value of the anisotropy. 

Now let us use our bosonization. The boundary action reads then, using the boundary conditions 
4'iL = (piR i^^ the UV to introduce the non chiral fields. 



J± 



/TT -^ I iTT ~ I IT ~ / iTT ~ 

2<?^i + Y Y</'2)(O)+rj;(TyCos(W-0i - \j—(t)2){Q) 



+ J//r^CT2Cos(V27r</.i(0)). (7.15) 



Consider now the expansion of the boundary free energy in powers of J^^Jjj. for every insertion of 
cos \/27r(^i, we need one insertion of cos( V^c^i + J^4'2) and one of cos(a/^(^i — 1/3^^2), contributing 
a term T^CTxTyCryTzCTz = —1 to the spin part. It follows that all th e spin terms actually disappear, and 
that the boundary free energy is the same as the one of the model ( [7.15 ) with no spin variables: 



Hint = -2J_icoaJ-(j)iCosJ — (i>2- J//cosV2n4>i. (7.16) 

This is nothing but the weak corrugation Hamiltonian of section 4 (eq. (4.17)), this time for a^ — ^, 
where the perturbation is marginal. 

Using the TBA for the anisotropic Kondo model, the impurity free energy is simply expressed as 



, d9 1 

-t^ irnr) — " -L 



■ln(l + e-'=^W/^), (7.17) 



™P" ^ ' 2^cosh(0-0B) 
where Tb oc e^ is the Kondo temperature, and ei is a pseudo energy, solution of the TBA system 

., ^ ,„ >,,e' -T^N,, I f JJ5^ m (1 + .-'">'-) ^ (7.S, 

Here, the labels k run over the nodes of the diagram 

1 2 4 t-3 / 

O — O (2) O — O t - 2 

t 



H.S. thanks P. Fendley and N. Warner for discussions on this point. 
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and we have restricted to the simplest values of the anisotropy such that -^^ = 1 — J// . 

Provided Jh > 0, the perturbation generates a flow (the situation here is quite different from the one 
boson theory, where the boundary cosine interaction is truly marginal) . Independently of the value of 
J||, one finds xi = e~'i("°°)/^ = 3, while yi = e-<°°'^/'^ = 2. The ratio ^ = -2=, the weU known 

value for the 4 channel Kondo model. This coincides with the ratio ^^ (4.30) for our values of the 
coupling constant. For Jn < (the equivalent of w^i < in section 4), meanwhile, the Kondo perturbation 
is irrelevant, and no flow is generated Hence, we conclude that for An — 1 (4.19), there is no new fixed 
point, and presumably the Y fixed point becomes identical with the N one in the limit a^ = ^. An 
e-expansion around this marginal point was developed by Kane and Fisher E 



VIII. CONCLUSION 

We have established the general connection between a simple model for QBM and boundary CFT 
with considerable attention paid to the issue of boson compactness. For a simple Hamiltonian we have 
conjectured a complete phase diagram, which is analogous to that of the boundary Potts model. These 
results were obtained by using a sort of "boundary embedding" involving the Potts model and the fusion 
approach to generating boundary states. We have also mentioned other fixed points which occur with 
more general Hamiltonians and which arise from fusion in the other sectors of the conformal embedding. 
We have discussed the integrability of some of the RG fiows. 

We note that the Y b.c. is more dynamical or " quan tum" than D or N. Within the weak corrugation 



formulation, while the potential term va of Eq. ( 4.17 ) is relevant, it does not simply "pin" the boson 



fields at one of its minima; rather they fiuctuate between the two inequivalent minima. This is only 
possible for a range of lattice spacing, 8tt/9 < a? <2tt where a semi-classical analysis fails. 

We have given an explicit demonstration that other, highly non-trivial b.c.'s are possible for 2 free 
periodic bosons besides D and N and their variants. However, since this demonstration rests on the 
conformal embedding and fusion, and since we don't have a conformal embedding for general values of 
the radius parameter, a, our construction only works at one special value of a. A more straightforward 
understandJHg of these b.c.'s, directly in terms of the bosons eludes us. Thus this general problem, known 
since 1992ffl remains open. 
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